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DEBRIS  FORMATION  AND  TRANSLATION  SUMMARY 


INTRODUCTION 

This  report  represents  a  systematic  effort  to  examine 
the  physical  basis  for  predicting  the  final  location  of  blast- 
initiated  debris.  There  are  three  principal  sources  of  this 
debris: 

•  Frangible  structural  elements,  such  as 
masonry  wall  panels. 

•  Nonfrangible  structural  elements,  such 
as  building  frames  and  wood  or  metal 
siding  and  roofing. 

•  Building  contents. 

The  first  two  categories  require  a  method  to  predict  the  loads 
at  which  they  will  come  apart  and  the  kinds  of  pieces  into  which 
they  will  break,  or  more  generally,  a  method  of  failure  predic¬ 
tion.  Chapters  Two  and  Three  of  this  report  deal  with  this 
problem. 

Presuming  a  knowledge  of  the  failure  modes,  the  impor¬ 
tant  question  from  a  postattack  point  of  view  is:  how  much  of 
these  elements  end  up  obstructing  the  adjacent  roadway?  More 
particularly,  there  is  interest  in  the  weight-size-composition, 
height,  and  total  volume  of  matter  in  the  desired  right-of-way. 
Chapter  Four  is  concerned  with  the  construction  of  a  computer- 
oriented  model  to  predict  the  distribution  of  "loose  particles", 
that  is,  structural  fragments  and  building  contents.  Also  in 
Chapter  Four,  assorted  loose  ends  are  tied  up  concerning  the 
finer  points  associated  with  the  transport  model. 

A  summary  of  the  state-of-the-art  in  debris  prediction 
is  shown  in  Table  1,  an  examination  of  which  will  show  that  with 
the  results  given  In  this  report,  the  theoretical  basis  for 
debris  prediction  is  pretty  well  covered.  However,  a  few  holes 
still  exist.  The  most  noticeable  is  the  restriction  of  the 
fragmentation  model  to  homogeneous  wall  panels.  Further  modi¬ 
fication  will  be  required  to  be  able  to  handle  nonhomogeneous 
wall  panels. 
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Since  a  great  number  of  walls  (including  those  made  of 


brick),  are  in  this  category,  such  an  extension  would  be  desir¬ 
able. 


Table  1 

SUMMARY  OF  DEBRIS- PREDICTION 


Debris  Source 

Method  of  Failure 
Prediction 

Method 

Location 

of  Final 
Prediction 

Frangible 

Fragmentation 

Transport  Model 

Structural  Elements 

Theory 

Nonfrangible 

Limited  Plasticity 

Continuity  (Frames) 

Structural  Elements 

Transport  Model 
(Siding  and  Roofing) 

Building  Contents 

Not  Applicable 

Transport  Model 
(Plus  overturning 
and  sliding  analysis 
for  diffraction- 
sensitive  items.) 

NONFRANGIBLE  STRUCTURAL  ELEMENTS 

Debris  resulting  from  the  effects  of  blast  on  nonfrang- 
ible  structural  elements,  such  as  beams  and  columns,  seems 
worthy  of  consideration  in  any  attempt  to  provide  meaningful 
inputs  for  postattack  recovery  planning.  This  follows  from  the 
fact  that  while  elements  of  this  sort  have  a  smaller  volume  of 
potential  debris  than  frangible  ones,  the  resulting  "particles" 
will  be  larger,  more  cumbersome,  and  hence,  more  demanding, 
pound  for  pound,  in  any  clean-up  effort.  With  this  motivation, 
we  have  striven  to  develop  an  analytical  procedure  capable  of 
predicting  the  size  and  weight  distribution  of  the  debris  de¬ 
posited,  in  a  nuclear  blast  environment,  by  elements  which  have 
some  ductility.  Such  elements  will  be  denoted  as  nonfrangible 
to  distinguish  them  from  frangible  (or  brittle)  ones,  such  as 
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unreinforced  wall  panels,  which  have  no  capacity  to  absorb 
energy  beyond  their  yield  points. 

For  all  practical  purposes,  the  load-response  behavior 
of  nonfrangible  structural  elements  can  be  divided  into  two 
categories,  based  on  the  plastic  regions  of  their  stress-strain 
diagrams.  The  response  is  either  sufficiently  ductile  to  allow 
the  use  of  an  elastic-perfectly  plastic  model  or  the  amount  of 
strain  that  can  be  accommodated  is  limited,  requiring  a  "limited 
plasticity"  model.  The  former  case,  which  has  been  thoroughly 
investigated  over  the  last  twenty  years,  is  generally  applicable 
to  steel- framed  structures.  The  latter  case,  which  is  appropri¬ 
ate  for  reinforced-concrete  structures,  was  considered  and  the 
effect  of  the  limited  ductility  was  demonstrated. 

Finally,  a  small  series  of  experiments  on  model  frames 
was  devised  to  check  the  validity  of  the  limited-plasticity  model 
and  verify  the  hypothesis  that  any  energy  supplied  to  a  frame  in 
excess  of  that  necessary  to  cause  collapse  is  taken  up  by  rota¬ 
tions  of  the  plastic  hinges  to  the  extent  of  their  capacities 
and  acceleration  of  the  mechanism,  rather  than  in  secondary  dam¬ 
age  between  hinges.  The  information  gained  from  this  series  of 
experiments  was  qualitative  in  nature. 

Some  conclusions  about  the  utility  of  the  theories  and 
techniques  demonstrated  are: 

•  The  limited-plasticity  theory  provides  a 
realistic  approach  for  predicting  blast- 
induced  debris  from  nonfrangible  structural 
elements  in  a  manner  which  Is  consistent 
with,  and  indeed  an  extension  of,  design 
procedures. 

•  Recourse  to  modern  computer-oriented  analysis 
techniques  overcomes  the  prohibitive  computa¬ 
tional  complexity  which  heretofore  has  inhib¬ 
ited  applications  of  limited  plasticity. 
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•  Models  of  reinforced-concrete  structures,  con¬ 
structed  at  low  cost  from  inexpensive  materials, 
can  be  used  to  provide  meaningful  answers  to 
questions  about  debris  production  which  charac¬ 
teristically  involve  gross  behavior  such  as 

the  collapse  mode. 

FRAGMENTATION  OF  FRANGIBLE  STRUCTURAL  ELEMENTS 

The  frangible  plate  structure  represents  a  significant 
debris  producing  element  in  the  form  of  wall  panels  and  a  vital 
source  of  dangerous  missiles  in  the  form  of  plate  glass.  The 
fragmentation  characteristics  of  such  structures  are  studied  in 
this  section  using  a  pragmatic  approach  which  blends  results 
from  statistical  fracture  theory  with  those  recently  obtained 
by  IITRI  on  an  experimental  study  of  dynamically- loaded  plaster 
plates  (Ref.  1).  The  work  extends  the  considerations  of  two 
previous  investigations  on  beam  fragmentation  to  the  plate 
(Ref.  2  and  3  ) . 

The  general  fragmentation  algorithm  consists  of  four 

steps: 

•  Determine  the  maximum  dynamic  stresses 
throughout  the  plate. 

•  Compute  the  probability  of  fracture 
initiation  throughout  the  plate. 

•  Divide  the  plate  into  appropriate  regions 
based  on  crack  propagation. 

•  Compute  the  distribution  of  fragment 
"sizes." 

Three  computational  procedures  are  described  for  deter¬ 
mining  the  distribution  of  fragment  sizes.  Each  of  these  methods 
begins  by  dividing  a  plate  into  regions  or  strips  formed  by  the 
principal  stress  trajectories.  These  strips  independently  frac¬ 
ture  or  remain  intact  and  the  combination  of  fracture  and 
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nonfracture  determines  the  geometry  and  number  of  fragments. 

The  first  computation  scheme,  the  combination  method,  considers 
individually  each  of  the  possible  2n  combinations  of  failure 
and  nonfailure  of  the  strips  where  n  is  the  total  number  of 
strips.  This  method  provides  the  specific  description  and 
quantity  of  every  possible  fragment,  and  in  addition,  it  details 
the  various  possible  mixtures  of  large  and  small  fragments.  It 
unfortunately,  is  very  time  consuming  even  with  the  aid  of  very 
large  computers. 

If  we  are  not  interested  in  how  the  various  fragments 
are  mixed  together,  we  can  adopt  a  very  efficient  procedure 
called  the  fragment  group  method,  for  calculating  the  total 
number  of  every  possible  type  of  fragment.  Here,  there  are 
only  (n/2)(n+i)  combinations  of  fragment  groups  to  be  consid¬ 
ered.  Although  the  increased  efficiency  of  the  method  of  frag¬ 
ment  groups  is  considerable,  an  even  faster  method  can  be  used 
if  we  again  settle  for  less  information.  The  final  method, 
calLed  the  method  of  runs,  determines  the  number  of  identical 
contiguous  nonfractured  strips.  It  will  not  furnish  information 
about  fragment  geometry;  only  fragment  weights. 

TRANSPORT  MODEL 

In  order  to  represent  the  effect  of  debris  transport  and 
subsequent  distribution,  it  is  necessary  to  move  from  a  problem 
space  consisting  of  the  real  world  to  a  more  abstract  mvthemati- 
cal  model.  This  abstraction  consists  of  representing  the  initial 
condition  of  possible  debris  as  a  series  of  lumped  masses  at  lev¬ 
els  above  ground.  Each  lumped  mass  is  characterized  by  a  unique 
particle  size  distribution.  The  particle  size,  in  turn,  has 
weight  and  shape  attributes  associated  with  it.  The  trajectory 
model  assumes  two  ideal  initial  conditions.  These  are: 

e  Zero  failure  time  of  fragmented  elements. 

e  An  initial  particle  velocity  of  zero. 

These  assumptions  were  made,  initially,  due  to  a  lack  of  know¬ 
ledge  concerning  any  other  possible  values. 
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A  study  concerning  these  parameters  has  since  been  made  and  is 
reported.  The  result  of  this  study  indicates  that  the  initial 
assumptions  are  well  grounded. 

SINBAD  (Simulation  Investigation  of  Nuclear  Blast  Asso¬ 
ciated  Debris)  is  a  problem- oriented  computer  language  that  deals 
with  the  problem  of  postattack  structural  debris.  In  a  previous 
investigation  (  Ref. 3  )  debris  profile  curves  (i,e.,  height  of 
debris  versus  distance  thrown)  were  developed  for  a  free-standing 
masonry  panel  wall.  Several  analyses,  both  manual  and  computer¬ 
ized,  were  utilized  to  predict  the  profile  of  a  single  wall.  The 
present  study  is  a  refinement  of  the  previous  techniques  and  is 
extended  to  include  any  grouping  of  walls  subjected  to  a  frontal 
shock.  It  is  now  also  possible  to  determine  the  size  distribu¬ 
tion  and  a  measure  of  the  momentum  of  the  debris  at  any  point  in 
the  profile.  The  language  is  expanoable  and  in  its  entirety  will 
include  frame  response  as  well  as  the  interior  contents  of  the 
structure.  The  flow  diagram  indicates  the  general  computational 
scheme.  The  boxes  that  are  now  dotted  are  components  that  will 
be  added  to  the  system  at  a  later  time. 
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Translation".  The  major  topics  investigated  were: 

•  Debris  generated  by  nonfrangible  structural 
elements . 

•  Fragmentation  of  plate- type  elements. 

•  Trajectory  of  debris  particles. 

Limited  investigations  were  also  performed  on  selected  topics 
relevant  to  debris  prediction. 
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ABSTRACT 

A  comprehensive  view  is  taken  of  the  physical  models 
required  to  estimate  volumes  and  heights  of  blast-initiated 
debris.  Particular  emphasis  and  development  is  directed  toward 
three  areas:  the  fragmentation  of  frangible  elements,  the  fail¬ 
ure  of  elements  with  limited  ductility,  and  the  transport  of 
debris  particles  by  blast  winds.  Computer  programs  to  handle 
the  computations  involved  in  these  three  models  have  been  written. 
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CHAPTER  ONE 
INTRODUCTION 


1.1  PERSPECTIVE  OF  DEBRIS  PREDICTION 

This  report  represents  a  systematic  effort  to  examine 
the  physical  basis  for  predicting  the  final  location  of  blast- 
initiated  debris.  There  are  three  principal  sources  of  this 
debris : 

•  Frangible  structural  elements,  such  as 
masonry  wall  panels. 

•  Nonfrangible  structural  elements,  such 
as  building  frames  and  wood  or  metal 
siding  and  roofing. 

•  Building  contents. 

The  first  two  categories  require  a  method  to  predict  the  loads 
at  which  they  will  come  apart  and  the  kinds  of  pieces  into  which 
they  will  break,  or  more  generally,  a  method  of  failure  predic¬ 
tion.  Chapters  Two  and  Three  of  this  report  deal  with  this 
problem. 

Presuming  a  knowledge  of  the  failure  modes,  the  impor¬ 
tant  question  from  a  postattack  point  of  view  is:  how  much  of 
these  elements  end  up  obstructing  the  adjacent  roadway?  More 
particularly,  there  is  interest  in  the  weight-size-composition, 
height,  and  total  volume  of  matter  in  the  desired  right-of-way. 
Chapter  Four  is  concerned  with  the  construction  of  a  computer- 
oriented  model  to  predict  the  distribution  of  "loose  particles", 
that  is,  structural  fragments  and  building  contents.  Also  in 
Chapter  Four,  assorted  loose  ends  are  tied  up  concerning  the 
finer  points  associated  with  the  transport  model. 

A  summary  of  the  state-of-the-art  in  debris  prediction 
is  shown  in  Table  1,  an  examination  of  which  will  show  that  with 
the  results  given  in  this  report,  the  theoretical  basis  for 
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Table  1 

SUMMARY  OF  DEBRIS -PREDICTION 


Method  of  Failure  Method  of  Final 
Debris  Source  Prediction  Location  Prediction 


Frangible 

Structural  Elements 

Nonfrangible 
Structural  Elements 


Building  Contents 


Fragmentation 

Theory 

Limited  Plasticity 
Not  Applicable 


Transport  Model 


Continuity  (Frames) 
Transport  Model 
(Siding  and  Roofing) 

Transport  Model 
(Plus  overturning 
and  sliding  analysis 
for  diffraction- 
sensitive  items.) 


debris  prediction  is  pretty  well  covered.  However,  a  few  holes 
still  exist.  The  most  noticeable  is  the  restriction  of  the 
fragmentation  model  to  homogeneous  wall  panels.  Further  modi¬ 
fication  will  be  required  to  be  able  to  handle  nonhomogeneous 
walls.  Since  there  are  a  great  number  of  walls  in  this  cate¬ 
gory,  which  Includes  those  made  of  brick,  such  an  extension* 
would  be  desirable. 
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CHAPTER  TWO 

DEBRIS  GENERATED  BY  NONFRANGIBLE  STRUCTURAL  ELEMENTS 


2.1  INTRODUCTION 

Debris  resulting  from  the  effects  of  blast  on  nonfrang- 
Ible  structural  elements,  such  as  beams  and  columns,  seems 
worthy  of  consideration  in  any  attempt  to  provide  meaningful 
inputs  for  postattack  recovery  planning.  This  follows  from  the 
fact  that  while  elements  of  this  sort  have  a  smaller  volume  of 
potential  debris  than  frangible  ones,  the  resulting  "particles" 
will  be  larger,  more  cumbersome,  and  hence,  more  demanding, 
pound  for  pound,  in  any  clean-up  effort.  With  this  motivation, 
we  have  striven  to  develop  an  analytical  procedure  capable  of 
predicting  the  size  and  weight  distribution  of  the  debris  de 
posited,  in  a  nuclear  blast  environment,  by  elements  which  have 
some  ductility.  Such  elements  will  be  denoted  as  nonfrangible 
to  distinguish  them  from  frangible  (or  brittle)  ones,  such  as 
unreinforced  wall  panels,  which  have  no  capacity  to  absorb 
energy  beyond  their  yield  points. 

For  all  practical  purposes,  the  load-response  behavior 
of  nonfrangible  structural  elements  can  be  divided  into  two 
categories,  based  on  the  plastic  regions  of  their  stress- strain 
diagrams.  This  distinction  is  shown  in  Fig.  1  for  a  bending 
member  where  moment  corresponds  to  stress  and  rotation  to  strain. 
The  response  is  either  sufficiently  ductile  to  allow  the  use  of 
an  elastic-perfectly  plastic  model  or  the  amount  of  strain  that 
can  be  accommodated  is  limited,  requiring  a  "limited-plasticity" 
model.  The  former  case,  although  rather  thoroughly  investigated 
over  the  last  20  years,  is  of  little  interest  for  debris-predic¬ 
tion  purposes.  The  latter  model,  however,  has  considerable  ap¬ 
plicability  (Ref.  1).  In  the  first  place,  the  removal  of  both 
metal  and  wooden  siding  from  building  frames  can  be  formulated 
as  a  limited  plasticity  problem  since  the  mode  of  failure  in¬ 
volves  both  rupture  at  connections  and  tearing  apart  of  the 
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Plastic  Moment 


Fig.  1  ELASTIC- PERFECTLY  PLASTIC  AND  LIMITED 
ROTATION  MODELS 


A 


panels  after  large  deformations.  This  aspect  was  not  treated 
in  detail  in  this  portion  of  the  research  effort,  but  rather, 
attention  was  focused  on  the  more  difficult  application  to  the 
collapse  of  reinforced  concrete  frames.  However,  the  goal  of 
this  study  was  an  ability  to  predict-  the  locations  on  a  rein¬ 
forced  concrete  frame  at  which  the  inelastic  rotations  will  be 
sufficiently  large  for  the  reinforcing  steel  to  be  exposed. 

The  steel  can  be  cut  at  these  points  and  the  structure  dismantled. 
This  corresponds  to  a  first  (or  ready-made)  level  of  debris 
clearance.  Attempts  to  break  the  structure  into  smaller  pieces 
will  require  chipping  away  the  concrete  in  order  to  sever  the 
reinforcement , 

2.2  LIMITED  PLASTICITY  OF  FRAMES 

The  behavior  of  reinforced  concrete  frames  is  still 
quite  a  controversial  subject  (Ref.  2).  One  point  cf  view  in¬ 
sists  that  for  multistory  frames  in  particular,  the  loss  of 
stiffness  due  to  the  beam-column  effect  must  be  considered  in 
ultimate  load  calculations  and  is  supported  by  experiments  on 
model  frameworks  (Ref.  3  through  6).  However,  an  analysis  of 
this  sort  neglects  the  support  given  by  walls  and  floors  and 
is  bound  to  be  overly  conservative  when  applied  to  complete 
buildings.  Even  in  an  examination  of  blast-load  effects  on 
framed  buildings,  where  the  walls  are  considered  to  have  been 
removed  by  the  diffraction  loading,  beam-column  effects  are 
important  only  in  tall,  slender  frames.  We  will  concentrate 
on  the  simpler  theory  which  is  applicable  to  the  great  majority 
of  buildings.  Another  basic  item  of  contention  is  the  choice 
of  a  model  to  represent  the  flexural  behavior  of  reinforced 
concrete.  One  side  demands  that  "strain-softening"  (i.e.,  re¬ 
sistance  increasing  to  a  maximum  and  then  decreasing  smoothly 
as  the  deformation  increases)  be  included  in  the  model.  (Consult 
the  papers  by  Barnard  and  Rosenblueth  in  Ref.  2).  The  other, 
and  preponderant,  viewpoint  expressed  in  the  paper  of  Baker  and 
Amarakone,  also  in  Ref.  2,  adopts  a  limited-plastic  model  of 

1 1 T  RESEARCH  INSTITUTE 


5 


V 


the  type  shown  in  Fig.  1.  We  have  gone  along  with  the  majority 
in  using  the  straight-line  model  for  limited  plasticity.  The 
reasoning  on  our  part  was  simple  since  we  are  concerned  with 
the  plastic  moment,  rotation  capacity,  and  energy  absorption 
(which  is  the  area  under  the  curve).  Thus,  whatever  the  exact 
shape  of  the  moment-rotation  curve  may  be,  a  la  Bernard  and 
Rosenbleuth,  we  can  pick  three  straight-line  segments  which  will 
match  those  salient  characteristics.  (This  is  done  at  the  ex¬ 
pense  of  accuracy  in  the  "elastic  rotation"  which  we  do  not  care 
about . ) 

Once  the  moment  at  a  section  becomes  equal  to  the  plas¬ 
tic  moment,  a  "plastic  hinge"  is  formed  and  the  rotation  in¬ 
creases  at  constant  moment.  If  we  postulate  a  limited  rotation 
capacity,  the  behavior  beyond  that  amount  is  like  a  "real  hinge" 
and  rotation  increases,  but  no  moment  is  transmitted  across  the 
section.  Clearly  then,  if  the  rotation  at  a  point  in  a  loaded 
structure  exceeds  the  capacity  at  that  point,  the  ultimate  load 
which  can  be  carried  will  be  the  same  as  would  be  indicated  by 
an  analysis  of  a  modified  structure  that  had  a  real  hinge  at  the 
point  in  question.  Moreover,  and  of  greater  interest  from  the 
debris  removal  aspect,  the  hinge  pattern  will  differ  in  general 
from  that  found  under  an  assumption  of  unlimited  rotation  capacity. 

In  order  to  have  the  ability  to  assess  rapidly  the  mag¬ 
nitudes  of  the  inelastic  rotations  encountered  in  a  large  framed 
structure,  a  computational  method,  first  suggested  by  Wang 
(Ref.  7)  as  a  limit  analysis  procedure,  was  programmed  for  IIT 
Research  Institute  (IITRI)  7094  computer.  Basically,  the  ap¬ 
proach  is  to  perform  a  sequence  of  elastic  analyses  on  the  struc¬ 
ture.  Consider  a  given  structure  and  loading  pattern.  If  an 
elastic  analysis  is  performed,  and  the  location  of  maximum  moment 
determined,  the  load  factor  can  be  adjusted  to  cause  a  plastic 
hinge  to  form  at  that  point.  Then,  after  adjusting  the  moments 
at  all  nodes  in  accordance  with  this  load  factor,  the  remaining 
moment  resistances  can  be  found.  Next,  an  elastic  analysis 
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performed  on  a  structure  tfhich  is  identical  with  the  original 
one  (except  for  a  pin  inserted  at  the  location  of  the  plastic 
hinge)  will  indicate  the  node  having  the  maximum  moment.  A 
load  factor  can  then  be  found  which  will  induce  a  moment  at  that 
point  equal  to  its  remaining  moment  resistance,  implying  the 
formation  of  a  plastic  hinge.  This  process  is  repeated  until  a 
collapse  mechanism  is  formed.  The  sum  of  the  load  factors  for 
all  cycles  is  the  ultimate  load  factor. 

This  method  may  seem  roundabout,  and  perhaps  it  is,  but 
it  is  well  suited  for  the  exceptionally  efficient  computer  solu¬ 
tions  utilizing  matrix  algebra,  since  the  modifications  can  be 
performed  automatically  during  the  analysis.  Also,  the  inelastic 
rotations  can  be  computed  at  each  stage,  permitting  inclusion 
of  the  effect  of  limited  rotation  capacity.  The  basis  for  the 
calculations  is  the  well-known  deflection  method,  where,  using 
matrix  notation,  the  member  end-rotations,  and  moments  £m], 

are  vectors  related  by  the  stiffness  matrix  S, 

(m)  -  i  (d)  (1) 

The  external  forces,  £  pj,  are  related  to  the  end  moments  by  the 
beam  and  bent  equations, 

[p]  >  a{m).  (2) 

It  follows  that  the  external  displacements,  /x],  are  related  to 
the  end  rotations  by 

^  i  }  -  AT{x},  where  AT  is  the  transpose  of  A.  (3) 

The  procedure  for  solution  is  to  use  Eqs.  (1),  (2),  and  (3)  to 
solve  for  the  displacements, 

[  X  }  -  I  A  S  AT]  1  [p],  (A) 
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rr, 


and  then  to  compute  the  end  moments  by 


[«} 


{x}. 


If  the  stiffness  matrix  of  the  original  structure  is 
designated  as  SQ,  the  end  rotations  at  the  end  of  the  first 
cycle  and  corresponding  to  the  formation  of  the  first  plastic 
hinge,  computed  by  either  Eq.  (1)  or  (3)  as 


[•>)  -  [Soj'U"} 


will  be  identical.  After  subsequent  cycles,  during  which  the 
stiffness  matrix  of  the  structure  has  been  modified,  the  in¬ 
elastic  rotation  at  the  nodes  H  will  be  given  by  the  difference 


(«}  -  [so]"1M  -aTW- 


A  listing  of  the  FORTRAN  IV  computer  program  is  given  in  Appen¬ 
dix  A. 

2.3  RESULTS  OF  FRAME  STUDIES 

A  limited-plasticity  analysis  of  a  framed  structure  car 
give  results  which  will  differ  from  those  of  a  standard  limit 
analysis  in  three  areas: 

•  the  ultimate  load  carried, 

•  the  total  energy  absorbed,  and 

•  the  final  collapse  mode. 

Two  example  problems  were  run  to  demonstrate  these  disparities 
The  first  example  demonstrates  the  reductions  in  both  ultimate 
load  and  energy-absorption  due  to  limited  rotation  capacity. 

The  frame  analysed  is  shown  in  Fig.  2  along  with  the  notation 
consistent  with  the  computer  program.  As  can  be  seen  in  Fig.  3, 
an  elastic-perfectly  plastic  (i.e.,  "limit")  analysis  will  in¬ 
dicate  collapse  at  a  load  factor  of  3000.  Now,  for  purposes  of 
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(b)  External  Force  - 
Displact'PU'nc 
Notation 


(c)  Internal  Moment  - 
Rotation 
Notation 


Fig.  2  FRAME  USED  IN  SAMPLE  PROBLEM 


I 


9 


members  have: 


Fig.  3  ELASTIC- PLASTIC  ANALYSIS  OF  ORIGINAL  FRAME 


illustration,  say  that  the  inelastic  rotation  capacity  at  all 
nodes  is  0.004  radian.  An  examination  of  the  computer  output 
for  this  problem,  which  is  displayed  as  Appendix  B,  will  show 
that  with  this  restriction  the  inelastic  rotation  at  node  8 
becomes  critical.  A  linear  interpolation  between  load  factors 
and  hinge  rotations  allows  us  to  fix  the  load  factor  consistent 
with  the  rotation  constraint  at  2630. 

To  find  the  energy  absorbing  capacity  of  the  frame,  we 
consider  the  load  which  can  be  supported  by  the  frame  under  im¬ 
posed  deformation.  When  the  rotation  capacity  at  node  8  is  ex¬ 
ceeded,  the  load  which  can  be  sustained  is  that  associated  with 
the  same  deformation  in  a  frame,  identical  with  the  original 
frame  except  for  a  real  hinge  at  node  8.  An  elastic -perfectly 
plastic  analysis  can  be  run  on  such  a  frame  and  the  effects  of 
the  rotation  constraints  found.  In  this  manner,  a  series  of 
modified  frames  can  be  considered  and  the  solid  curve  shown  in 
Fig.  4  constructed.  The  area  under  this  curve  is  a  measure  of 
the  energy  which  can  be  absorbed  by  the  frame  in  question. 

To  illustrate  the  possibility  of  restrictions  on  rotation 
capacity  leading  to  a  different  collapse  mode,  the  frame  shown 
in  Fig.  5,  6,  and  7  was  analyzed.  Since  node  16  proved  to  be 
critical  in  this  case,  it  was  assumed  that  its  rotation  capacity 
would  be  exceeded  while  that  of  all  other  nodes  would  not.  The 
structure  was  then  analyzed  with  a  real  hinge  inserted  at  node  16. 
The  results  are  shown  in  Fig.  8  and  9.  It  can  be  seen  that  not 
only  is  the  collapse  load  lowered,  but  also  the  mode  of  collapse 
differs,  since  dead  loads  are  included.  In  the  previous  example, 
since  only  side-on  loads  acted,  the  collapse  had  to  be  in  a  side 
sway  mode. 
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action  In  Direction 


Fig.  6  POSITIVE  DIRECTIONS  FOR  END 
MOMENTS  AND  ROTATIONS 


Fig.  7  POSITIVE  DIRECTIONS  FOR  EXTERNAL 
FORCES  AND  DISP.  ACEMENTS 


COLLAPSE  MECHANISM  FOR  CASE  NO.  2, 
p  *  840 


2.4  QUALITATIVE  FRAME  EXPERIMENTS 

A  small  series  of  experiments  on  model  frames  was  devised 
t  ’  check  the  validity  of  the  limited-plasticity  model  and  verify 
the  hypothesis  that  any  energy  supplied  to  a  frame  in  excess  of 
that  necessary  to  cause  collapse  is  taken  up  by  rotations  of  the 
plastic  hinges  to  the  extent  of  their  capacities  and  acceleration 
of  the  mechanism  rather  than  in  secondary  damage  between  hinges 
The  information  to  be  gained  from  this  series  of  experiments  was 
qualitative  in  nature. 

The  geometry  of  the  frames  tested  is  shown  in  Fig.  10. 

The  materials  used  were  Hydrostone  plaster  and  a  soft  wire  rein¬ 
forcement.  The  mold  used  to  case  the  frames  is  displayed  in 
Fig.  11.  Due  to  the  small  percentage  of  reinforcement,  about 
1  percent,  the  behavior  of  the  frames  was  governed  almost  entirely 
by  the  reinforcement.  Static  collapse  load  predictions  are  shown 
in  Fig.  12  and  the  observed  collapse  loads  in  Fig.  13  and  14. 

Since  the  objectives  were  qualitative  in  nature,  the  static  col¬ 
lapse  tests  were  performed  in  a  Riehle  testing  machine  for  ease 
of  load  application.  The  fact  that  the  load  scale  on  this  ma¬ 
chine  only  permitted  readings  to  the  nearest  10  lb  was  still 
sufficient  to  show  satisfactory  agreement  between  prediction 
and  observation.  Further  verification  of  the  limited-plasticity 
theory  was  gained  from  the  static  collapse  test  on  the  single¬ 
story  frame.  The  history  of  the  failure  was  as  follows: 

•  At  a  load  somewhat  below  100  lb,  cracks 
became  visible  at  the  column  bases. 

•  Deformation  continued  without  increase  in 
load  at  about  100  lb. 

•  As  deformation  increased,  the  load  fell 
suddenly  to  about  50  lb. 

•  After  further  deformation  at  this  level, 
the  load  fell  to  zero. 
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Since  the  moment  due  to  the  compressive  force 
is  negligible. 


Fig.  12  STATIC  COLLAPSE  LOAD  PREDICTION 
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Fig.  14  STATIC  COLLAPSE,  TWO- STORY  FRAME 
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This  sequence  of  events  is  consistent  with  a  limited-plasticity 
formulation  shown  in  Fig.  15. 

Finally,  dynamic  collapse  tests  were  performed  on  both 
single-story  and  two-story  frames.  Since  the  behavior  of  frames 
in  a  high-yield  blast  environment  is,  for  all  intents  and  pur¬ 
poses,  solely  dependent  on  response  to  drag  loadings  of  durations 
much  greater  than  the  natural  period  of  the  structure,  a  dynamic 
loading  fixture  was  devised  to  produce  a  load  pulse  as  shown  in 
Fig.  16.  Loads,  both  slightly  greater  than  the  observed  static 
collapse  loads  and  more  than  twice  as  much,  were  applied  in  this 
fashion.  The  collapse  modes  and  amounts  of  damage  at  the  hinges 
were  comparable  in  all  cases,  as  was  predicted.  The  responses 
of  the  four  frames  tested  under  impact  are  shown  in  Fig.  17 
through  20. 

2.5  CONCLUSIONS 

Some  conclusions  about  the  utility  of  the  theories  and 
techniques  demonstrated  in  this  chapter  are  appropriate: 

•  The  limited-plasticity  theory  provides  a 
realistic  approach  for  predicting  blast- 
induced  debris  from  nonfrangible  structural 
elements  in  a  manner  which  is  consistent 
with,  and  indeed  an  extension  of,  design 
procedures . 

•  Recourse  to  modern  computer-oriented  analysis 
techniques  overcomes  the  prohibitive  computa¬ 
tional  complexity  which  heretofore  has  inhib¬ 
ited  applications  of  limited  plasticity. 

•  Models  of  reinforced-concrete  structures,  con¬ 
structed  at  low  cost  from  inexpensive  materials, 
can  be  used  to  provide  meaningful  answers  to 
questions  about  debris  production  which  charac¬ 
teristically  involve  gross  behavior  such  as 
the  collapse  mode. 
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Fig.  15  LIMITED -PLASTICITY  BEHAVIOR  OF  TEST  FRAME 


Side-on 


Deflection  at 
complete  static 
collapse 


Fi*.  16  DYNAMIC  LOADIN'^  OF  TEST  FRAMES 


C  COLLAPSE,  SINGLE -STORY  FRAME 


Fig.  20  DYNAMIC  COLLAPSE,  TWO- STORY  FRAME 
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CHAPTER  THREE 
PLATE  FRAGMENTATION 


3.1  INTRODUCTION 

The  frangible  plate  structure  represents  a  significant 
debris  producing  element  in  the  form  of  wall  panels  and  a  vital 
source  of  dangerous  missiles  in  the  form  of  plate  glass.  The 
fragmentation  characteristics  of  such  structures  are  studied  in 
this  section  using  a  pragmatic  approach  which  blends  results 
from  statistical  fracture  theory  with  those  recently  obtained 
by  IITRI  on  an  experimental  study  of  dynamically- loaded  plaster 
plates  (Ref.  8) .  The  work  we  shall  describe  extends  the  con¬ 
siderations  of  two  previous  programs  on  beam  fragmentation  to 
the  plate  (Ref.  9  and  10) . 

In  the  first  of  these  programs,  the  statistical  nature 
of  the  problem  is  established  together  with  the  physical  as¬ 
sumptions  underlying  the  basic  computational  scheme.  Essentially, 
the  method  considers  separately  every  possible  combination  of 
crack  patterns.  As  such,  it  provides  a  description  of  the  dis¬ 
tribution  of  fragment  shapes  and  masses,  and  in  addition,  it 
can  be  used  to  characterize  the  mixture  of  different  fragments. 
Unfortunately,  the  computational  time  for  this  program  is  very 
great  even  for  large  computers.  In  the  second  beam  fragmenta¬ 
tion  program,  a  very  efficient  and  rapid  computation  method 
related  to  the  theory  of  runs  was  proposed  which  described  only 
the  fragment  size  distribution  -  the  original  locations  of  the 
fragments  cannot  be  determined  nor  are  they  required  for  beam 
response.  As  we  shall  see,  this  additional  information  may  be 
useful  for  describing  the  fragmentation  of  plates. 

The  general  fragmentation  algorithm  consists  of  four 

steps : 

•  Determine  the  maximum  dynamic  stresses 
throughout  the  plate. 
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•  Compute  the  probability  of  fracture 
initiation  throughout  the  plate. 

•  Divide  the  plate  into  appropriate  regions 
based  on  crack  propagation. 

•  Compute  the  distribution  of  fragment 
"sizes." 

Each  of  these  steps  is  discussed  in  the  following  subsections. 

3.2  DYNAMIC  STRESS  ANALYSIS 

To  decide  whether  or  not  fracture  will  initiate  at  a 
point  in  a  dynamically  loaded  plate,  we  must  first  know  the 
"worst"  stress  state  that  can  occur  at  the  point.  This  is  a 
straightforward  determination  when  no  fractures  occur  throughout 
the  load  history.  If,  on  the  other  hand,  fractures  do  develop 
during  the  loading  process,  the  problem  is  considerably  more  com¬ 
plicated.  Even  for  a  material  with  a  deterministic  strength  we 
would  have  to  consider  changing  boundary  conditions,  the  speed 
of  crack  growth,  and  the  direction  of  crack  propagation.  For  a 
brittle  material  with  statistically  distributed  strength,  the 
number  of  combinations  requiring  analysis  would  truly  be  enormous. 

To  extricate  ourselves  from  this  forbidding  prospect, 
we  have  introduced  the  assumption  that  the  maximum  dynamic 
stresses  are  independent  of  the  fracture  characteristics  of  the 
structure.  The  following  comments  are  relevant  to  this  approxi- 
mat ion : 

1.  No  experimental  evidence  has  been  sought  to 
examine  the  validity  of  this  assumption  for 
different  types  of  dynamic  loading. 

2.  The  unloading  that  accompanies  the  first 
fracture  oi  a  slowly  loaded  statically 
determinate  beam  usually  precludes  a 
second  fracture. 
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3.  Multiple  fractures  invariably  occur  on  a 
rapidly  loaded  statically  determinate  beam. 

4.  The  more  severe  the  dynamic  loading  the 
smaller  the  fragment  size  and  the  greater 
the  number  of  fragments. 

5.  Under  such  an  assumption  the  various  possible 
fracture  patterns  are  stochastically  independent. 

6.  Crack  velocity  is  substantially  below  the  veloc¬ 
ity  of  elastic  disturbances. 

7.  The  actual  stress  magnitudes  in  a  structure  will 
usually  be  equal  to  or  lower  than  those  computed 
for  a  dynamically  equivalent  plate  with  infinite 
strength.  This  implies  that  we  will  experience 
fewer  crack  initiations  and  larger  pieces  than 
we  might  predict. 

Consistent  with  our  principal  assumption  of  independence, 
i.e.  5  (above),  we  shall  proceed  to  calculate  the  maximum  stresses 
occurring  in  a  rectangular  simply- supported  plate  subject  to  uni¬ 
form  load  across  its  surface  but  varying  in  time.  The  coordinate 
system  and  plate  dimensions  are  shown  in  Fig.  21.  Conventional 
small  deformation  theory  is  used  and  the  plate  is  assumed  to  be 
homogeneous  and  isotropic . 

As  a  specific  example,  we  have  chosen  a  simply- supported 
rectangular  plate  with  an  exponentially  decaying  load, 
q  38  pQ  exp  pt].  The  initial  velocities  and  displacements 
are  taken  to  be  zero.  The  deflection  for  such  a  plate  is  de¬ 
scribed  in  Section  9.5  of  Ref.  11  where  their  general  deflection 
expression,  Eq.  (8),  can  be  specialized  by  taking  q  *  pQexp(-|j,t) 
and  f  *  g  *  0.  Then,  using  Eq.  (10)  of  this  reference,  we  ob¬ 
tain  after  a  simple  integration  the  required  plate  deflection: 
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12(l-v2)7 


h  *  plate  thickness 
E  *  Young's  modulus 
v  *  Poisson's  Ratio 
7  *  weight  density 

g  -  acceleration  of  gravity. 

The  resulting  moments  can  be  found  by  substituting  Eq.  (8) 
into  the  following  which  relate  moments  to  deflections. 

.3 
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-Eh' 


xx 
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-Eh' 


yy 
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12(l-v  ) 
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r  2  2 
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Lay  J 
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dxdy 

(c) 
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12(l-v") 

It  is  then  possible  to  find  the  principal  moments  from: 

MPM2  *  1/2  Kcx  +  Myy]  ■ 


(10) 
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Since  the  principal  stresses  are  related  to  the  principal 
moments  by 

Si  2  »  “j  M^  2  »  (11) 

we  can  find  the  magnitudes  and  directions  of  the  principal 
stresses  in  the  plate  at  any  time. 

Because  of  the  arduous  summations  involved,  Eqs.  (8) 
through  (11)  were  programmed  for  the  IBM  7094  digital  computer. 

A  particular  problem,  that  of  a  square  plate  15  in.  on  a  side 
and  1/2  in.  thick,  was  run  and  the  resulting  contours  of  maximum 

principal  stress  are  shown  in  Fig.  22  for  *  5  psi  and 

-1  ° 

U  =  2  sec  .  The  curves  are  the  contours  at  the  time  when  the 

stress  at  the  center  of  the  plate  (which,  of  course,  is  the  max¬ 
imum  stress  in  a  simply-supported  plate)  is  a  maximum,  i.e., 
t  -  0.001958  sec.  The  maximum  stresses  are  very  closely  approx¬ 
imated  by  the  stresses  associated  with  the  contour  lines  in 
Fig.  22  because  the  plate  deflects  predominantly  in  the  first 
mode. 

3.3  PROBABILITY  OF  FRACTURE  INITIATION 

Using  the  principal  stresses  calculated  by  the  methods 
of  the  previous  subsection,  we  shall  address  ourselves  to  the 
problem  of  establishing  the  probability  that  fracture  will  ini¬ 
tiate  in  a  typical  subdivision  of  the  plate  shown  in  Fig.  23. 
These  subdivisions  are  identified  by  the  integers  running  from 
1  to  120  and  their  associated  bending  moments  are  calculated  at 
their  centroids.  Figure  24  shows  a  subdivision  from  which  we 
have  extracted  a  slice  which  is  subjected  to  the  principal 
stresses  (SpS2»0).  Before  we  can  establish  its  reliability, 
it  is  necessary  that  a  theory  be  developed  for  multiaxial  stress 
fields. 
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i  «  b  -  15  Inches 
i  ■  1/2  Inch 


Fig.  22  LINES  OF  CONSTANT  MAXIMUM  STRESS 


3.3.1  Combined  Stress  Theory 


In  his  classic  paper  of  1939  (Ref.  12),  Weibull  developed 
an  expression  for  the  fracture  probability  of  a  brittle  material 
under  a  polyaxial  stress  state.  Using  a  different  point  of  view, 
we  shall  expand  on  his  brief  statistical  treatment  of  this  com¬ 
bined  stress  problem  and  extend  our  results  to  cases  of  varying 
mechanical  and  thermal  loading,  and  to  materials  which  cannot 
be  represented  by  the  Weibull  distribution  function. 

Briefly,  it  is  our  objective  to  establish  a  fracture 
surface,  i.e.,  to  find  a  relationship  among  the  strengths  achieved 
under  various  stress  states.  The  usual  approach  to  this  problem 
in  either  brittle  or  ductile  materials  is  to  find  a  property 
common  to  all  stress  states  that  will  indicate  failure  or  non¬ 
failure.  In  ductile  materials  the  distortion  energy  represents 
such  a  property,  since  incipient  flow  occurs  in  any  stress  state 
in  which  the  distortion  energy  is  equal  to  the  distortion  energy 
obtained  in  a  tension  specimen  at  yield.  Stated  in  another  way, 
we  can  correlate  yielding  under  any  stress  state  with  the  dis¬ 
tortion  energy.  Our  approach  for  brittle  materials  is  completely 
analogous -we  shall  try  to  find  a  property  that  will  correlate 
with  the  reliabilities  associated  with  the  various  possible  com¬ 
bined  stress  conditions. 

To  avoid  the  "size  effect"  problem  observed  in  the 
strength  of  brittle  elements,  (i.e.,  increasing  fracture  stress 
with  decreasing  volume)  we  shall  begin  our  study  by  considering 
a  finite  unit  volume  AV  of  fixed  size.  We  assume  that  both 
the  material  and  the  stress  state  in  this  unit  volume  are  homo¬ 
geneous  and  that  the  materials  used  in  all  the  unit  volumes  to 
be  considered  have  been  drawn  from  the  same  population.  In  ad¬ 
dition,  we  shall  restrict  the  study  to  brittle  materials  that 
are  statistically  isotropic,  i.e.,  the  distribution  of  strengths 
obtained  from  an  indefinitely  large  number  of  unit  volumes  will 
be  identical  in  every  direction. 
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Fig.  23  PLATE  SUBDIVISIONS  SHOWING  THEIR  RISK  OF  RUPTURE  VALUES 


24  TYPICAL  PLATE  SUBDIVISION 


We  shall  assume  that  the  principal  stresses  S^,  S2, 
which  act  on  a  basic  unit  volume  are  proportional  to  a  load 
factor  S,  that  is, 

■  aS 

S2  -  @S  (12) 


where  a,  £,7  are  constants  which  define  the  stress  state.  The 
strength  of  a  basic  element  will  be  taken  as  the  maximum  load 
factor  that  it  can  equilibrate.  Failure  of  the  unit  element  is 
represented  by  its  inability  to  equilibrate  the  applied  loading. 
It  is  important  to  point  out  that  it  is  possible  for  cracks  to 
initiate  and  propagate  within  the  unit  volume  without  causing 
failure  of  the  element.  Materials  in  which  cracks  can  be 
arrested  or  which  provide  alternative  load  paths  when  local 
failures  occur  are  classified  as  parallel  or  series-parallel 
materials.  If  a  local  failure  necessarily  leads  to  overall 
failure,  the  associated  material  is  called  a  series  or  '‘weakest 
link"  material.  One  can  advantageously  adopt  an  infinitesimal 
unit  volume  for  the  series  material  and,  as  we  shall  subse¬ 
quently  discuss,  combined  stress  testing  is  greatly  simplified 
in  this  case. 

Only  the  tensile  or  cohesive  mode  of  failure  will  be 
considered  in  this  investigation.  We  shall  assume  that  neither 
compressive  nor  shear  stresses  influence  the  strength  of  a 
brittle  material.  The  potential  usefulness  of  this  tension 
criterion  is  a  consequence  of  two  observations;  first,  that 
the  shear  strength  of  brittle  materials  is  usually  an  order  of 
magnitude  greater  than  the  tensile  strength,  and,  second,  that 
it  is  extremely  difficult  to  eliminate  tensile  stresses  from 
prototype  or  laboratory  elements.  Almost  every  structural 
failure  of  a  brittle  component  can  be  attributed  to  the  pres¬ 
ence  of  some  distribution  of  tensile  stresses. 
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3.3.2  Two-Dimensional  Theory  Heuristic  Development 


When  we  attempt  to  describe  the  statistical  fracture 
strength  of  a  unit  volume  of  material  under  a  uniaxial  stress 
state,  the  axial  stress  (strain)  is  the  only  reasonable  choice 
for  the  statistical  variate.  Taking  a  general  form  for  any  cumu¬ 
lation  distribution  function,  we  can  write  the  fracture  probabil¬ 
ity  F  for  the  uniaxial  stress  state  as 


F(a)  -  1  -  exp 


g(a) 


where  AV  is  the  specified  volume  of  the  basic  unit  element,  v 
is  a  volume  of  unity,  and  o  is  the  axial  stress.  The  delinea¬ 
tion  of  the  constant  AV/v  does  not  affect  the  generality  of 
this  expression  and  in  the  special  case  of  a  series  material  it 
provides  a  convenient  representation.  If  we  examine  the  strength 
of  a  unit  volume  of  an  isotropic  material  under  a  general  homo¬ 
geneous  stress  state,  it  follows  that  failure  will  depend  only 
on  the  three  principal  stresses  acting  on  the  unit.  Thus,  the 
probability  of  failure  of  the  unit  volume  can  be  designated  as 
F(Si, 82,83)  where  the  three  principal  stresses  are  taken  as  the 
statistical  variates.  For  this  case  we  shall  take  Eq.  (13)  in 
the  form 


-<n  (1-f<si.s2.s3>] 


g ( s 1 , S2 , S3) 


(14) 


For  a  specified  reliability  (1-F),  we  note  that  Eq.  (14)  becomes 
g(S-pS2»S2)  “  constant,  which  defines  our  fracture  surface. 

On  the  basis  that  failure  is  caused  only  by  tensile 
stresses,  it  seems  reasonable  to  look  for  the  function  g  with¬ 
in  the  collection  of  all  possible  tensile  stresses  which  can  oc¬ 
cur  at  any  point  in  the  unit  volume.  In  the  plane  stress  prob¬ 
lem  we  can  relate  the  normal  stress  an  acting  in  any  direction 
to  the  principal  stresses  through  the  expression 


cos^0  +  S2  sin^'0 


(15) 
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where  0  is  the  angle  between  an  and  S^.  As  0  sweeps  through 
all  values  from  -tt/2  to  7t/2,  Eq.  (15)  describes  every  possible 
normal  stress  acting  at  a  point.  The  normal  stresses  associated 
with  the  various  directions  described  by  0  are  shown  in  Fig.  25 
for  several  different  stress  states.  The  question,  now,  is 
what  are  the  distinguishing  features  of  these  figures  that  will 
reflect  the  differences  they  cause  in  a  material's  response? 

The  most  obvious  first  guess  is  to  differentiate  among 
these  stress  states  by  comparing  the  areas  associated  with  the 
tensile  normal  stresses.  However,  this  approach  does  not  re¬ 
flect  the  possibility  that  the  magnitude  of  the  stresses  may 
have  a  different  influence  than  their  extent  or  distribution. 

For  example,  hydrostatic  and  pure  tension  stress  states  are 
depicted  in  Fig.  26  that  lead  to  the  same  area  but  where  one 
peak  stress  is  twice  the  other.  Experience  indicates  that  the 
pure  tension  state  is  the  more  critical.  This  suggests  that  we 
"weight"  the  ordinates  in  these  figures  and  then  compare  the 
areas  of  the  weighted  normal  stress- theta  diagrams. 

Assuming  a  statistically  isotropic  material,  the  weight¬ 
ing  should  be  independent  of  the  orientation  0  of  the  normal 
stress.  We  might  use  for  example  a  power  function  to  modify  the 

normal  stresses,  i.e..  Do  where  D  and  k  are  constants.  This 

n 

alteration  results  in  the  dashed  curve  shown  on  the  left  side 
of  Fig.  27.  If  the  normal  stress  distribution  for  several 
stress  states  were  weighted  in  this  fashion,  we  could  compare 
the  areas  of  the  resulting  curves,  that  is, 

itSj.Sj)  •  Aim  .  D  f  oj  *  (16) 

n~ 

where  the  integration  extends  over  those  values  of  0  where  the 
normal  stress  is  non-negative.  Because  of  symmetry  we  need  con¬ 
sider  only  the  positive  normal  stresses  in  the  interval  zero  to 
7i/2.  To  account  for  the  possibility  that  tensile  stresses  below 
a  certain  magnitude  a ^  may  not  cause  failure,  we  may  choose  to 
weight  the  difference  (CTn"^p  as  shown  in  the  right  half  of 
Fig.  27. 
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Pure  Tension 


Hydrostatic  Tension 


Biaxial  Tension 


Pure  Shear 


Fig.  25  NORMAL  STRESS  DISTRIBUTIONS  FOR  VARIOUS  STRESS  STATES 


The  associated  area  is  given  by 


g(SvS2)  -  Area  -  D  J  ( on-c£)k  d0  - 


an  »ai 


Certainly,  the  use  of  a  power  function  to  weight  the 
normal  stress- theta  diagrams  is  completely  arbitrary  and  there 
are  many  other  ways  of  manipulating  and  distorting  such  curves. 
Our  problem  is  to  find  a  weighting  function  that  will  reflect 
the  influence  of  stress  state  on  the  reliability  of  a  unit 
volume.  Denoting  the  weighting  function  by  f,  the  fracture 
probability  becomes 

F(SltS2)  -  1  -  exp  -  — §¥  J  f(vV  d9  ]•  (18) 

°n*°l 


We  are  now  in  a  position  to  describe  certain  guidelines 
for  the  selection  of  f.  First,  to  account  for  the  possible 
existence  of  a  zero  fracture  probability  stress  Oj ,  we  must  take 


f<v°P 


f  •  0 


°n—  at-° 


The  latter  condition  implies  that  both  and  an<* 

that  in  such  cases  F  ■  0.  At  the  other  extreme,  we  expect  that 
fracture  is  a  certainty  when  either  or  S2  is  positive  and 
unbounded;  hence,  F  •  1  implies  that 

f  00  when  Sj-*  +  «. 
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Furthermore,  we  would  expect  on  physical  grounds  that  the 
failure  probability  would  increase  continuously  with  increasing 
principal  stresses,  thus, 


f  ...  continuous  and  monotone  increasing. 


Finally,  f  must  be  chosen  in  such  a  way  that  the  associated 
F(SpS2)  fits  the  cumulative  distribution  curve  obtained  from 
fracture  tests  conducted  using  various  stress  states.  In 
particular,  it  is  necessary  that  fracture  data  obtained  under 
pure  tension  be  represented  by  F(S^,0). 

Typical  examples  of  admissible  forms  for  the  weighting 
function  f  are  the  following: 


f  « 


(19) 


f  - 


exp  [a(cn-aj 


1 


(20) 


f 


(21) 


f  -  +  b (<v°|)2  ♦  c<V°i)3  +  ••• 

A> 0,  B20,  C»0 


(22) 


where  a,  k,  A,  B,  C,  aQt  and  Oj  are  constants  of  the  material. 
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3.3.3  Three-Dimensional  Theory 

The  extension  of  our  theory  given  in  Eq.  (18)  to  three 
dimensions  requires  that  we  appropriately  distort  the  surface 
formed  by  the  normal  stress  vector  in  three  dimensions.  This 
vector  is  given  in  the  polar  coordinates  as 

on  -  cos^Sjeos2*  +  S2sin2^)  +  S3sin2<(  (23) 

where  the  angles  4  and  ^  are  defined  in  Fig.  28a.  A  typical 
surface  representing  the  focus  of  normal  stress  vectors  is 
shown  in  Fig.  28b  for  a  biaxial  tension  field.  A  weighted 
surface  is  formed  by  and  its  volume  can  be  introduced 

into  the  general  distribution  function  to  give 


F(S^,S2»S3)  ■  1-exp 


[-  ?  jlf  “V 

L  ar£?l 


oB)dV  . 


(24) 


Specializing  to  the  form  of  f  given  in  Eq.  (19)  and  using  polar 
coordinates ,  the  failure  probability  F  is  given  by 

r  2  rrfn  „  _  3k  n 

F  -  1-exp  -  -j  ~  d*  cosrf  drf  (-A ,J,)  (25) 


where  we  find  three  cases: 
1)  si  i  S2  >  Qt 


*  cos’ 


SjCOs  iM^sin  * 
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2)  S1  >  ar  S2  <  ajL 


3)  si  <  <*j>  s2  f 

*  0  (F  *  °)  • 

Equation  (25)  can  be  written  in  the  form 

1  -  F  *  e'B  (26) 

where  the  "risk  of  rupture"  B  is  given  by  the  negative  of  the 
term  within  the  square  brackets  of  Eq.  (25).  The  risk  of  rupture 
B  was  evaluated  numerically  for  each  slice  of  every  plate  sub¬ 
division  indicated  in  Fig.  23.  Specifically,  the  following 
data  was  used: 

Plate  size:  15  x  15  x  1/2  in. 

Overpressure:  P0  =  5  psi 

Pressure  decay:  \i  =  2  sec 

Statistical  parameters:  k  ■  3 

a0  *  1500  psi 
o^  *  50  psi 

Now,  a  value  of  B^  for  the  ith  slice  shown  in  Fig.  24 
enables  us  through  Eq.  (26)  to  establish  the  probability  that  no 
fracture  will  initiate  in  the  slice,  (1-F^).  The  probability 
that  no  fracture  will  initiate  in  the  entire  subdivision  (1-Fg) , 
requires  the  simultaneous  survival  of  each  slice,  thus, 

n 

(1-FS)  -  (1-F1)(l-F2)...(l-Fn)  -  TT  <1‘Fi)  (27) 

;  m\  b 
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where  n  is  the  total  number  of  slices.  Substituting  Eq.  (26) 
into  this  equation  we  obtain 

(1-FS)  =  exp(-Bg)  =  exp  (-  ^  Bt)  .  (28) 

Therefore,  the  risk  of  rupture  of  a  "big  piece"  is  equal  to  the 
sum  of  the  risk  of  ruptures  of  its  component  "small  pieces". 

The  sum  of  the  slice  risk  of  ruptures  for  each  plate  subdivision 
is  tablulated  in  Table  2,  together  with  its  centroid  coordinates 
and  maximum  principal  bending  stresses.  These  risk  of  rupture 
values  are  displayed  in  Fig.  23  by  the  lower  number  in  each 
subdivision. 

3.4  PLATE  EXPERIMENTS 

One  of  the  most  difficult  ispects  of  the  plate  fragmenta¬ 
tion  problem  concerns  the  question  of  crack  propagation.  Crack 
initiation  was  the  concern  of  the  previous  two  subsections.  In 
the  beam  problem,  when  a  crack  initiated  within  the  beam  volume 
this  always  resulted  in  a  fracture  surface  which  was  roughly 
perpendicular  to  the  beam  axis.  When  a  crack  initiates  within  a 
plate,  its  direction  of  travel  is  not  obvious.  Furthermore,  we 
meet  a  new  problem  when  many  cracks  are  propagating  because 
one  crack  crossing  the  path  of  a  second  crack  will  generally 
arrest  the  second  crack.  We  are  faced,  therefore,  with  the  "who 
got  there  first"  problem.  In  the  face  of  these  complications, 
we  examined  the  results  of  experiments  conducted  with  Hydrostone 
plaster  plates  under  dynamic  loadings.  The  experiments  con¬ 
ducted  at  IITRI  (Ref.  8 ) ,  were  supposed  to  demonstrate 
characteristic  crack  patterns  that  would  provide  the  needed 
propagation  information  for  our  fragmentation  analysis.  If  no 
patterns  were  obtained,  our  analysis  procedure  would  have  to  be 
abandoned,  and  indeed,  the  hope  of  developing  a  rational  prediction 
scheme  would  be  pretty  gloomy.  Fortunately,  patterns  did  emerge 
from  these  tests  and  we  shall  very  briefly  summarize  the  findings 
which  are  described  in  detail  in  Ref.  8. 
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3.4.1  Description  of  Drop  Test 


It  has  been  shown  in  Ref.  8  that  the  response  of  a  plate 
under  any  uniform  time -dependent  loading  can  be  made  identical 
to  that  achieved  in  a  drop  test  when  the  appropriate  support  de¬ 
celeration  is  imposed.  To  produce  the  dynamic  load  in  our  drop 
test  facility,  the  plate  support  was  mounted  on  the  drop  table 
as  shown  in  Fig.  29.  The  idea  was  to  drop  the  table  and  suddenly 
decelerate  it,  which  would  load  the  plate  mounted  on  the  supports 
(as  shown  in  Fig.  30)  with  downward  acting  inertia  body  forces. 

To  increase  the  downward  loading,  sand  was  piled  onto  the  plate 
and  held  in  place  by  the  box  device  in  Fig.  31.  The  results  of 
a  typical  drop  test  are  illustrated  in  Fig.  32  where  the  frag¬ 
ments  are  held  intact  by  masking  tape  on  their  upper  surface. 

To  check  out  the  symmetry  of  the  drop  test  leading,  two 
plastic  plates  were  stress  coated  and  dropped  from  different 
heights.  As  can  be  observed  from  Fig.  33  and  34,  the  loading 
is  excellent  and  a  pattern  of  principal  directions  is  obtained 
which  is  not  unlike  that  obtained  for  the  pressure  loading 
q  *  pc  exp  [-ut]  as  shown  in  Fig.  22. 

3.4.2  Results  and  Conclusions 

Typical  examples  of  the  crack  outlines  obtained  for  five 
different  size  Hydrostone  plaster  plates  are  shewn  in  Fig.  35. 

We  first  observe  that  these  cracks  form  a  pattern.  Second,  by 
comparing  the  crack  pattern  on  the  square  plates  to  the  stress 
coat  patterns  of  Fig.  33  and  34  we  see  that  for  the  most  part 
the  cracks  propagate  along  the  principal  direction  trajectories. 
Further  examination  of  the  square  plates  indicates  that  the  cen¬ 
tral  pattern  forms  first.  In  all  of  the  cases,  cracks  occur 
along  45  deg  lines  at  the  corners. 

On  the  basis  of  these  observations,  we  shall  postulate 
the  formation  of  the  primary  fracture  mode  shown  in  Fig.  3Ca  and 
the  secondary  fracture  mode  shown  in  Fig.  36c.  The  strips  in 
the  secondary  mode  are  intended  to  approximate  the  principal  stress 
trajectories.  MT  RESearch  institute 
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Table  2 

PRINCIPAL  STRESSES  AND  RISKS  OF  RUPTURE 
(P  ”  5  psi;  u  «  2  sec-1) 
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Figure  35  provides  typical  fracture  patterns  of  rectan¬ 
gular  panels  with  length-to-width  ratios  of  approximately  2,  3, 
and  4.  The  fracture  patterns  are  generally  what  would  be  ex¬ 
pected.  As  the  length  to  width  ratio  of  the  plate  increases, 
the  performance  of  the  plate  appears  to  approach  that  of  one 
supported  on  the  two  long  sides  only.  The  "square"  center  sec¬ 
tion  of  a  square  plate  associated  with  the  primary  failure  mode 
apparently  rather  rapidly  degenerates  from  a  square  through  a 
rectangular  phase  and  into  essentially  a  line.  Figure  35,  for 
example,  tends  to  indicate  that  for  even  a  length-to-width  ratio 
of  2,  the  center  section  has  almost  entirely  degenerated.  Thus, 
the  prediction  of  the  primary  fracture  mode  for  rectangular 
plates  may  be  simpler  than  for  square  plates.  It  would  appear 
to  follow  from  the  degeneration  of  the  plate’s  center  region  to 
a  line  that  debris  fragment  sizes  might  be  derived  on  the  basis 
of  the  procedures  for  the  secondary  fracture  mode  alone. 
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Fig.  31  PLASTER  PLATE  TEST  WITH  SAND  OVERBURDEN 
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Fig.  32  TYPICAL  PLATE  FRAGMENTATION 


3.5  FRAGMENTATION  ANALYSIS 


Before  we  describe  the  methods  for  defining  primary  and 
secondary  fragments  wj  must  first  establish  the  probability  that 
a  crack  will  initiate  within  an  arbitrary  region  of  a  plate.  Con¬ 
sider  the  triangular  area  enclosed  by  OBC  for  the  square  plate 
of  Fig.  23.  Assume  that  the  principal  stress  trajectories  sug¬ 
gest  that  the  area  be  subdivided  into  the  four  strips  identified 
by  the  Roman  numerals  shown  in  the  figure.  Note  that  strip  II 
contains  the  plate  subdivisions,  8,  9,  10,  11,  12,  24,  25,  26, 

27,  37,  38,  39,  40,  41,  51,  52,  53,  54,  62,  63,  64  and  65.  The 
risk  of  rupture  for  strip  II,  is  equal  to  the  sum  of  the 

risks  of  rupture  associated  with  the  preceding  sequence  of  sub¬ 
divisions,  i.e., 

BH  *  Bg  +  Bg  +  ...  +  B64  +  B65  *  0.00558 

Now,  the  probability  that  a  crack  will  initiate  in  strip  II  is 
simply 

Fu  =  -  i.,-0.00558 

We  may  now  consider  the  primary  mode. 

3.5.1  Primary  Fracture  Mode 

To  establish  the  size  of  the  central  square  fracture 
pattern  we  will  divide  the  central  region  of  the  plate  into  the 
imaginary  square  strips  shown  in  Fig.  36b.  The  failure  proba¬ 
bility  F^  of  each  of  the  strips  will  be  computed  and  the  length 
associated  with  the  largest  F^  or  B^^  will  be  taken  as  the 
size  of  the  square  pattern. 

In  Fig.  37  we  have  computed  one- eighth  of  the  risk  of 
rupture  for  each  of  the  square  strips  shown  in  Fig.  36b.  The 
maximum  occurs  in  the  strip  containing  the  subdivisions  111,  112, 
113,  and  114.  Clearly  then,  this  defines  the  primary  fracture 
mode.  We  note  in  passing  that  the  maximum  stresses  decrease  as 
we  move  away  from  the  center  and  that  the  primary  strip  volumes 
increase  as  we  move  from  the  center.  This  explains  why  we  find 
a  relative  maximum  between  the  center  and  the  edges. 


Fig.  37  RISKS  OF  RUPTURE  ALONG  HORIZONTAL  STRIPS,  PRIMARY  MODE 


Ig.  37  RISKS  OF  RUPTURE  ALONG  HORIZONTAL  STRIPS,  PRIMARY  MODE 
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We  observed  from  the  various  replications  of  drop  tests 
on  square  plates  that  the  size  of  the  center  square  remained 
fairly  constant.  If  large  variations  would  have  occurred  the 
probability  of  getting  a  size  is  simply  F^,. 

3.5.2  Secondary  Fracture  Mode 

Using  the  hypothesis  that  the  secondary  cracks  will  fol¬ 
low  the  principal  stress  trajectories,  we  can  divide  the  trape¬ 
zoidal  regions  formed  by  the  primary  cracks,  into  the  fracture 
strips  shown  in  Fig.  38.  Each  of  these  strips  will  independently 
fracture  or  remain  intact  in  exactly  the  same  manner  previously 
described  for  beam  fragmentation.  Specifically,  the  simply  sup¬ 
ported  plate  gives  rise  to  the  same  problem  solved  for  the  simply 
supported  beam. 

There  are  three  methods  available  for  dealing  with  secon¬ 
dary  plate  fractures.  We  begin  each  method  by  numbering  the 
strips  as  shown  in  Fig.  38.  Then,  using  the  same  technique  em¬ 
ployed  to  find  the  probability  of  fracture  initiation  in  strip  II 
of  Fig.  23,  we  can  find  the  fracture  probability  of  each  of  the 
n  strips  indicated  in  Fig.  36,  Fm.  Finally,  we  observe  that  the 
periphery  ABCD  is  a  free  boundary.  We  may  now  consider  each 
method  separately. 

3 . 5 . 2 . 1  Combination  Method 

This  method,  which  is  described  in  Ref.  9,  considers  in¬ 
dividually  each  of  the  possible  2n  combinations  of  failure  and 
nonfailure  of  the  strips.  If  F^  is  the  fracture  probability  of 
the  ith  strip  and  Sj^  the  associated  svirvival  probability  (note: 
S|  ■  1-F^),  then  the  following  combinations  of  fracture  and  sur¬ 
vival  tabulated  in  Table  3  are  possible  in  a  four-strip  plate. 
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Table  3 

POSSIBLE  COMBINATIONS  IN  A  FOUR-STRIP  TRAPEZOID 


Combinations 


S1S2S3S4 

F1S2S3S4 

S1F2S3S4 

S1S2F3S4 


S1S2S3F4 

F1F2S3S4 

F1S2F3S4 

f1S2S3F4 


S1F2F3S4 

S1F2S3F4 

S1S3F3F4 

F1F2F3S4 


FXF2S3F4 

F1S2F3F4 

SXF2F3F4 

F1F2F3F4 


Each  of  these  products  represent  the  probability  that  the  repre¬ 
sented  combination  will  occur.  The  sum  of  these  probabilities 
will,  of  course,  equal  unity. 

Now,  let  us  examine  a  typical  combination,  say  the  under¬ 
lined  one,  and  describe  its  significance  to  the  fragmentation 
problem.  First,  if  n  plates  are  dynamically  loaded,  4n  trapezoids 
will  give  rise  to  secondary  fractures.  Consequently,  the  number 
of  times  the  underlined  combination  will  occur  is  AntSjFjS^S^) . 
Associated  with  this  particular  combination  is  the  mixture  of  the 
twc  fragments  (strip  1)  and  (strips  3  +  4).  An  examination  of 
Table  3  indicates  that  these  two  fragments  can  arise  from  other 
combinations;  for  example,  strip  (3  +  4)  is  formed  by  both 

S1F2S3F4  anf*  F1F2S3S4*  is  *  8*nPl®  matter  of  bookkeeping  to 
accumulate  the  number  of  times  each  possible  fragment  occurs.  On 
the  other  hand,  it  is  very  time  consuming  to  consider  each  of  the 
2n  possible  combinations  which  generate  the  various  fragments. 

The  type  and  efficiency  of  debris  removal  equipment  will 
be  influenced  in  a  significant  way  by  the  composition  of  the  de¬ 
bris.  By  studying  the  more  frequently  occurring  combinations  of 
fracture  and  nonfracture,  it  is  possible  to  estimate  the  charac¬ 
ter  of  a  mixture  of  fragments.  The  combination  which  appears 
most  frequently  is  associated  with  the  following  probability. 
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(29) 


pm»x*  TT 

If  4  0.5  this  combination  is  unique. 


3. 5. 2. 2  Fragment  Group  Method 


If  we  are  not  Interested  in  how  the  various  fragments 
are  mixed  together,  we  can  adopt  a  very  efficient  procedure  for 
calculating  the  total  number  of  every  possible  type  of  fragment. 
In  an  n-strip  plate  segment  there  are(n/2) (n+l)possible  combina¬ 
tions  of  contiguous  strips.  Each  of  these  combinations  repre¬ 
sent  a  possible  fragment.  We  can  easily  display  these  combina¬ 
tions  as  shown  in  Table  4  for  a  four-strip  trapezoid.  The 
fragments  are  designated  by  the  numbers  of  the  strips  contained 
in  the  fragment.  For  example,  fragment  2,3  is  composed  of  the 
strips  2  and  3  in  Fig.  38.  To  obtain  this  fragment,  it  is  clear 
from  this  figure  that  strips  1  and  4  must  fracture,  and  strips  2 
and  3  must  not  fracture.  The  probability  of  this  happening  is 
represented  as  the  probability  of  simultaneously  getting  fracture 
in  strips  1  and  4  and  getting  no  fracture  in  strips  2  and  3,  i.e., 


FlS2S3F4* 


As  another  example,  we  see  that  fragment  1,2  can  be  rea¬ 
lised  by  survival  of  strips  1  and  2  followed  by  fracture  in 
strip  3.  It  does  not  matter  whether  strip  4  fractures  or  not. 
Thus,  the  probability  of  obtaining  fragment  1,2  in  a  trapezoid 
is  simply  (S^F^  *  l).  The  total  mnsber  of  fragments  "1,2" 
realised  from  n  plate  experiments  is  4n(S^$2F2>. 


If  a  fragment  is  composed  of  strips  k,  k+1,  ...,k+/,  the 
probability  of  its  occurrence  in  a  trapezoid,  Pk  ^  , 

is  given  by 


i*k+< 


Fk,k+l, . . .  ,k+*  “  Fk-1  Fk+*+l  fl*Fi> 


(30) 


where  the  fracture  probabilities  F0  and  Fn+i  represent  the  imagi¬ 
nary  strips  shown  in  Fig.  38. 
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For  a  trapezoid  with  a  free  boundary,  Fo  -  F„+r  l-  • 
fixed  boundary  condition  is  represented  by  FQ“  1.  Until  the 
fracture  pattern  is  known,  we  cannot  comment  on  the  shape,  loca¬ 
tion  or  behavior  of  the  n+1  strip  in  the  fixed  boundary  plate. 

Table  4 

NUMBER  AND  TYPE  OF  FRAGMENTS  IN  FOUR- STRIP  TRAPEZOID 


Fragment 

Designation 

Probability  of 
Occurrence-Trapezoid 

1 

S1F2 

2 

W3 

3 

F2S3F4 

4 

F3S4 

1,2 

S1S2F3 

2,3 

F1S2S3F4 

3,4 

F2S3S4 

1,2,3 

S1S2S3S4 

2,3.4 

FlS2S3S4 

1.2, 3,4 

SlS2S3S4 

3, 5.2. 3  Method  of  Runs 

By  considering  every  one  of  the  possible  2n  distinct 
fracture  patterns,  the  method  of  combinations  provides  the  specif 
ic  description  and  quantity  of  every  possible  fragment,  and  in 
addition,  it  details  the  various  possible  mixtures  of  large  and 
small  fragments.  The  method  of  fragment  groups  sacrifices  this 
latter  Information,  but  it  Increases  the  computational  efficiency 
enormously.  For  example,  if  the  number  of  strips  n  is  equal  to 
20,  the  combination  method  considers  2n  •  1,048,576  distinct 
fracture  patterns;  the  method  of  fragment  groups  will  consider  at 
most  all  of  the  possible  fragment  combinations,  (n/2)(n+l)*2010. 


-m 


Although  the  increased  efficiency  of  the  method  of  frag¬ 
ment  groups  is  considerable,  an  even  faster  method  can  be  used 
if  we  settle  for  less  information.  This  method,  called  the  meth¬ 
od  of  runs,  was  described  in  Ref.  7  for  the  fragmentation  of 
beams.  The  procedure  as  developed  is  not  directly  applicable  to 
the  plate  problem.  To  see  this  we  shall  consider  the  general 
problem  of  describing  the  fragments  resulting  from  fractures  in 
strips  1  and  3.  For  the  beam  we  would  say  that  we  had  a  one-unit 
piece  between  units  1  and  3  and  between  unit  3  and  the  support; 
that  is,  units  2  and  4  remain  intact.  Therefore,  for  this  com¬ 
bination  we  would  have  recorded  two  "one-unit"  pieces.  In  the 
plate,  a  glance  at  Fig.  38  indicates  that  strips  2  and  4  are 
different  and  we  cannot  claim  generally  that  we  have  two  one- 
unit  strips. 

With  this  in  mind,  we  shall  begin  our  approach  by  selec¬ 
ting  strips  with  equal  areas  as  shown  in  general  in  Fig.  39a  or 
in  particular  in  Fig.  39b.  Now,  every  one-strip,  two-strip,  and 
r-strip  fragment  (or  run)  has  the  same  mass. 

This  choice  of  equal  strip  areas  reduces  the  fragmenta¬ 
tion  problem  exactly  to  that  described  for  beams  in  Ref. 10.  For 
example,  to  find  the  total  number  of  two-unit  runs  in  the  four 
strips  shown  in  Fig.  39b,  we  observe  that  a  two-unit  strip  can 
occur  in  two  ways:  two  nonfailure  strips  followed  by  a  failure 
at  either  end  or  two  nonfailure  preceded  and  followed  by  fail¬ 
ures.  The  probability  that  these  events  will  take  place  is 
given  by: 

(1-Fj)(1-F2>Fj. 1  two-unit  fragment,  bottom 

(l-F^Hl-Fj^.  1  two-unit  fragment,  top 

two-unit  fragment,  middle 

The  sum  of  these  individual  probabilities  is  the  total  probabil¬ 
ity  of  obtaining  a  two-unit  fragment  from  one  trapesoid  in  the 
plate. 
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Fig.  39  GEOMETRIC  PROPERTIES  FOR  EQUAL  STRIP  AREAS 


The  general  formula  for  calculating  the  probability  that 
a  run  of  r  equal  area  strips  will  occur  in  a  trapezoid  is  given 
by 

n-£ 

P(r)  =  Z  Pk>k+1, . . .  ,k+£  (31> 

k*l 

or 

n-r+1  k+r-1 

P(r)'I  Fk-lFk+r  TT  <32> 

k*l  i=k 

where  FQ  and  F  are  the  fracture  probabilities  in  the  two  imag¬ 
inary  strips  shown  in  Fig.  38.  Here,  F0*Fn+^=  1.  Computing  the 
fragmentation  from  this  formula  is  very  rapid  and  inexpensive; 
however,  we  know  only  the  weight  characteristics  of  the  fragments, 
not  their  geometry  or  their  mixture. 

As  a  final  comment  we  should  note  that  the  propagation 
cf  a  crack  is  at  best  a  temperamental  and  sensitive  phenomenon. 

One  should  not  be  surprised  if  a  single  crack  branches  into  two 
cracks,  or  if  an  occasional  crack  propagates  across  the  princi¬ 
pal  stress  trajectories.  These  pecularities  will  produce  a  larger 
number  of  small  fragments  and  a  smaller  number  of  large  fragments 
than  predicted. 
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CHAPTER  FOUR 

TRAJECTORY  OF  DEBRIS  PARTICLES 

4.1  DESCRIPTION  OF  THE  PHYSICAL  MODEL 

In  order  to  represent  the  effect  of  debris  transport  and 
subsequent  distribution,  it  is  necessary  to  move  from  a  problem 
space  consisting  of  the  real  world  to  a  more  abstract  mathemati¬ 
cal  model.  This  abstraction  consists  of  representing  the  initial 
condition  of  possible  debris  as  a  series  of  lumped  masses  at  lev¬ 
els  above  ground.  Each  lumped  mass  is  characterized  by  a  unique 
particle  size  distribution.  The  particle  size,  in  turn,  has 
weight  and  shape  attributes  associated  with  it.  The  trajectory 
model  assumes  two  ideal  initial  conditions.  These  are: 

•  Zero  failure  time  of  fragmented  elements. 

•  An  initial  particle  velocity  of  zero. 

These  assumptions  were  made,  initially,  due  to  a  lack  of  know¬ 
ledge  concerning  any  other  possible  values.  A  study  concerning 
these  parameters  has  since  been  made  and  is  reported  at  the  con¬ 
clusion  of  this  chapter.  The  result  of  this  study  indicates  that 
the  initial  assumptions  are  well  grounded. 

4.2  INTRODUCTION  TO  SINBAD 

SINBAD  (Stimulation  investigation  of  Nuclear  Blast  Asso¬ 
ciated  Debris)  is  a  problem-oriented  computer  language  that  deals 
with  the  problem  of  postattack  structural  debris.  In  a  previous 
investigation  (Ref. 10)  debris  profile  curves  (i.e,,  height  of 
debris  versus  distance  thrown)  were  developed  for  a  free-standing 
masonry  panel  wall.  Several  analyses,  both  manual  and  computer¬ 
ized,  were  utilized  to  predict  the  profile  of  a  single  wall.  The 
present  study  is  a  refinement  of  the  previous  techniques  and  is 
extended  to  include  any  grouping  of  walls  subjected  to  a  frontal 
shock.  It  is  now  also  possible  to  determine  the  size  distribu¬ 
tion  and  a  measure  of  the  momentum  of  the  debris  at  any  point  in 
the  profile.  The  language  is  expandable  and  in  its  entirety  will 
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include  frame  response  as  well  as  the  interior  contents  of  the 
structure.  The  flow  diagram  illustrated  in  Fig.  40  indicates 
the  general  computational  scheme.  The  boxes  that  are  now  dotted 
are  components  that  will  be  added  to  the  system  at  a  later  time. 
The  remaining  sections  of  this  chapter  describe  the  input  lan¬ 
guage  and  sample  problems  run  on  the  program. 

4.2.1  Input  Language 

The  form  of  the  input  to  the  SINBAD  processor  differs 
significantly  from  most  other  computer  programs.  Format  and 
ordering  of  card  input  have  been  almost  eliminated;  they  have 
been  replaced  by  a  set  of  commands  consistent  with  postattack 
terminology.  The  fact  that  a  group  of  characters  starts  with  a 
letter  is  sufficient  to  recognize  a  word.  Similarly,  a  number 
indicates  numerical  data;  a  decimal  point  distinguishes  a  deci¬ 
mal  number  from  an  integer;  and  a  blank  or  a  comma  after  a  group 
of  characters  indicates  the  end  of  the  group. 

The  input  commands  may  be  data  descriptors,  data  to  be 
stored,  or  more  generally  information  about  the  input  process. 

A  data  descriptor  (e.g.,  YIELD  or  OVERPRESSURE)  communicates  to 
the  system  that  the  number  that  follows  is  to  be  associated  with 
that  command.  Data  to  be  stored  consist  of  the  numerical  data 
associated  with  data  descriptors.  Commands  such  as  WEAPON  PARAM¬ 
ETERS,  PREBLAST  STRUCTURAL  CONFIGURATION  and  SOLVE  actually  con¬ 
trol  the  internal  flow  of  the  program.  Table  5  contains  the 
dictionary  of  available  commands.  Each  command  occupies  a  sep¬ 
arate  input  card  in  the  data  and  a  card  may  be  continued  by  plac¬ 
ing  a  dollar  sign  ($)  in  the  first  column  of  the  following  cards. 
Each  input  card  is  printed  on  the  system  output  before  the  solu¬ 
tion  phase  of  the  processor  takes  over.  It  is  possible  to  put 
comment  cards  into  the  input  phase  simply  by  placing  an  asterisk 
(*)  in  column  1  of  the  card.  This  card  is  simply  echo  printed, 
but  otherwise  ignored.  Table  6  illustrates  a  set  of  commands 
that  is  sufficient  to  describe  a  debris  problem.  Once  the  prob¬ 
lem  has  been  initially  described  for  one  wall  and  solved,  it  is 
necessary  to  change  only  those  parameters  which  one  wishes  to 
vary  in  any  subsequent  wall  or  problem. 
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Table  5 

DICTIONARY  OF  PROCESS  COMMANDS  AND  DATA  DESCRIPTORS 


Process  Commands 
WEAPON  PARAMETERS 

PREBLAST  STRUCTURAL 
CONFIGURATION 


FRAGMENTATION  CHARACTERISTICS 


COMPUTE  FRAGMENTATION 
CHARACTERISTICS 


OUTPUT 


SOLVE 


Data  Descriptor 

YIELD 

OVERPRESSURE 
GROUND  ZERO  DISTANCE 

WALL  HEIGHT 
HEIGHT  BETWEEN  FLOORS 
SPACE  BETWEEN  WALLS 
NORMALIZING  FACTOR 

NUMBER  OF  PARTICLE  SIZES 
PARTICLE  SIZES 
PERCENTAGE  BY  SIZE 
ACCELERATION  COEFFICIENT 

BEAM  REPRESENTATION 

LENGTH 

DEPTH 

WIDTH 

STRESSO 

STRESSU 

PROFILE  DISTRIBUTION 
LOCATIONS 

DISTANCES  FROM  FIRST  WALL 
VELOCITY  DESCRIPTION 
DEBRIS  PROFILE  PLOT 


Table  6 

A  SUFFICIENT  SET  OF  COMMAS  AND  INPUT  TO 
SPECIFY  A  DEBRIS  PROBLEM  TO  SINBAD 


WEAPON  PARAMETERS 

YIELD  5000  KT 

OVERPRESSURE  10  PSI 

PREBLAST  STRUCTURAL  CONFIGURATION 

WALL  HEIGHT  40  FLOORS 

HEIGHT  BETWEEN  FLOORS  10  FEET 

DISTANCE  OF  WALL  FROM  INITIAL  WALL  50  FEET 

NORMALIZING  FACTOR  1.0 

FRAGMENTATION  CHARACTERISTICS 

NUMBER  OF  PARTICLE  SIZES  5 

PARTICLE  SIZE?  2.0,  4.0,  6.0,  6.0,  10.0  INCHES 

PERCENTAGE  BY  SIZE  0.05,  0.32,  0.16,  0.32,  0.05 

ACCELERATION  COEFFICIENT  0.0 

OUTPUT 

PROFILE  DISTRIBUTION 
DISTRIBUTION  OF  SIZES 
LOCATIONS  3 

DISTANCES  FROM  FIRST  WALL  50,  150  ,  300  FEET 
DEBRIS  PROFILE  PLOT 
SIZE  DISTRIBUTION  PLOT 
SOLVE 

WEAPON  PARAMETERS 
OVERPRESSURE  -  20.0  PSI 
SOLVE 


The  command  SOLVE  terminates  the  input  phase  of  the  processor 
and  transfers  control  to  the  computational  section.  When  the 
specified  problem  is  solved  and  the  answer  printed,  control  is 
automatically  returned  to  the  input  phase.  Each  of  the  data 
descriptors  will  now  be  discussed  in  detail. 

e  The  process  command  WEAPON  PARAMETERS  has  three 
data  descriptors:  YIELD,  OVERPRESSURE,  and  GROUND 
ZERO  DISTANCE.  The  YIELD  is  the  weapon  size  in 
kilotons,  and  is  used  in  conjuction  with  either 
the  OVERPRESSURE  (psi)  or  GROUND  ZERO  DISTANCE  (ft) 
to  specify  an  overpressure-distance  relationship. 

This  relationship  is  presently  based  on  a  mach 
region  surface  burst  assumption,  however,  as  the 
overall  system  is  modular  in  concept,  airburst  and 
regular  reflection  capabilities  could  be  included 
with  only  some  additional  effort.  Again,  it  is 
only  necessary  to  specify  either  the  OVERPRESSURE 
or  the  GROUND  ZERO  DISTANCE.  Knowledge  of  one  of 
these  parameters,  along  with  YIELD,  is  sufficient 
to  determine  the  other. 

e  PREBLAST  STRUCTURAL  CONFIGURATION  consists  of  four 
data  descriptors  that  describe  the  wall  under  in¬ 
vestigation.  WALL  HEIGHT  gives  the  total  number 
of  floors  (i.e.,  panels)  in  the  wall.  HEIGHT  BETWEEN 
FLOORS  is  the  panel  height  in  feet.  SPACE  BETWEEN 
WALLS  is  the  distance  in  feet  of  the  wall  presently 
being  investigated  from  the  last  previously  investi¬ 
gated  wall.  If  only  one  wall,  or  the  initial  wall 
in  a  multiwall  configuration  is  being  investigated, 
this  descriptor  is  unnecessary.  Finally,  a  NORMALIZ¬ 
ING  FACTOR  descriptor  is  Included  to  account  for  the 
normalization  of  the  debris  profile  curve. 
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This  normalization  has  been  explained  in  the  previ¬ 
ous  report  (Ref. 10)  and  it  suffices  to  say  that  this 
descriptor  is  usually  the  product  of  an  individual 
panel's  length  and  thickness.  If  the  NORMALIZING 
FACTOR  is  unity,  then  the  subsequent  debris  profile 
will  be  normalized  by  a  unit  width  volume  (i.e.,  the 
product  of  the  length  and  thickness,  sq  ft,  of  an 
individual  panel) . 

The  process  commands  FRAGMENTATION  CHARACTERISTICS 
of  COMPUTE  FRAGMENTATION  CHARACTERISTICS  describe 
the  type  of  particles  that  result  due  to  panel 
fragmentation.  This  report  will  only  include  a 
description  of  the  FRAGMENTATION  CHARACTERISTICS 
process  command  since  the  computational  model  of 
panel  fragmentation  is  only  in  a  formative  stage 
at  present.  The  panel  fragmentation  model,  based 
on  a  beam  analogy  that  was  developed  in  the  previous 
report,  has  een  included  in  the  present  system  but 
has  net  been  utilized.  This  was  done  because  its  use 
was  considered  marginal  in  light  of  the  work  done 
on  panel  fragmentation  as  discussed  in  Chapter  Three. 
Thus,  the  data  descriptors  listed  under  COMPUTE 
FRAGMENTATION  CHARACTERISTICS  are  consistent  with  the 
input  necessary  for  that  previous  fragmentation  anal¬ 
ysis.  The  data  descriptor  NUMBER  OF  PARTICLE  SIZES 
indicates  the  number  of  different  size  particles  re¬ 
sulting  from  panel  fragmentation.  PARTICLE  SIZES  is 
the  descriptor  of  an  array  of  the  individual  particle 
sizes  in  inch  units  and  each  is  separated  by  a  comma. 
This  array  is  listed  in  descending  order  of  size. 

In  a  similar  manner,  PERCENTAGE  BY  SIZE  is  a  corres¬ 
ponding  array  of  the  percentages  of  an  individual 
panel  falling  into  each  of  the  previously  described 
particle  sizes.  ACCELERATION  COEFFICIENT  describes 
the  shape  and  orientation  in  flight  of  an  individual 
debris  projectile. 


Under  normal  usage  this  parameter  is  set  equal  to 
zero  and  the  program  assumes  a  sphere  of  an  equiv¬ 
alent  volume  radius.  If  it  is  desired  to  investi¬ 
gate  other  shapes  with  several  orientations,  then 
ACCELERATION  is  equal  to  2  x  mas 8 /projected  area 
in  units  of  lb/sq  in. 

•  The  process  descriptor  OUTPUT  controls  the  type 
of  printed  and  computed  output  that  can  be  obtained 
from  the  system.  PROFILE  DESCRIPTION  indicates 
that  a  record  of  debris  height  as  a  function  of 
distance  from  an  initial  wall  is  desired.  The 
next  two  data  descriptors  are  utilized  to  obtain 
the  percent  by  size  range  at  each  desired  location 
in  the  debris  profile.  LOCATIONS  is  the  number  of 
points  in  the  debris  profile  where  a  size  distribu¬ 
tion  breakdown  is  wanted.  DISTANCES  FROM  FIRST  WALL 
is  the  array  of  distances  in  feet  from  the  initial 
wall  to  the  points  in  the  debris  profile  where  a 
size  distribution  is  desired.  VELOCITY  DESCRIPTION 
generates  three  output  relationships:  cumulative 
debris  momentum,  minimum  debris  momentum,  and  maxi¬ 
mum  debris  momentum  as  a  function  of  distance  from 
the  initial  wall  in  feet.  These  relationships  are 
normalized  by  the  mass  of  an  individual  panel  and 
are  actually  momentum  per  unit  length  along  the 
debris  profile.  This  will  be  further  developed  in 
the  following  section.  The  data  descriptor  DEBRIS 
PROFILE  PLOT  results  In  the  machine  plotting  of 
the  different  relations  previously  outputed  in 
printed  form. 

e  The  process  command  SOLVE  transfers  control  from 
the  input  phase  to  the  computational  mode. 

The  process  commands  and  data  descriptors  described 
above  and  in  Table  5  can  be  lnputed  In  any  order,  however,  the 
order  outlined  in  Table  6  seems  to  '«  logical. 
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Once  a  problem  has  been  initially  defined  it  is  only  necessary  co 
re-spec ify  those  data  descriptors  that  change  in  any  subsequent 
problem.  This  is  also  illustrated  in  Table  6  by  the  change  in 
overpressure  from  10  to  20  psi. 

4.3  MOMENTUM  ANALYSIS 

One  of  the  primary  effects  of  nuclear  associated  struc¬ 
tural  debris  is  the  tertiary  effect  it  has  on  the  unsheltered 
population.  It  has  been  shown  (Ref.  13)  that  whereas  one  may 
survive  from  free- field  prompt  effects  of  a  nuclear  explosion, 
(i.e.,  blast,  thermal  and  radiation)  he  may  still  be  highly  vul¬ 
nerable  to  high-speed  flying  debris  projectiles.  In  order  to 
measure  the  effectiveness  of  this  type  of  phenomenology  the  pro¬ 
jectile's  mass  as  well  as  its  speed  must  be  included.  This  is 
accomplished  by  describing  the  projectile's  momentum  per  unit 
length  over  the  debris  profile. 

The  trajectory  analysis  that  was  utilized  to  find  the 
final  position  of  flying  debris  also  yields  the  projectile's 
final  speed.  Figure  41  illustrates  how  a  normalized  momentum 
per  unit  length  is  determined.  A  ire  range  is  specified  by  two 
projectile  sizes.  Each  of  these  sizes  has  a  final  speed  associ¬ 
ated  with  It  as  well  as  a  final  horizontal  displacement  from  its 
original  position.  The  normalized  momentum  per  unit  length  is 
determined  by: 


I  V, 

t*i  - 


(33) 


where 


is  the  norms lized  debris  momentum  per 
unit  length 

is  the  percent  of  debris  falling  into 
size  range  1  which  is  composed  of 
particle  sizes  falling  between  sizes 
J  and  j+l 
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Pig.  41  FORMULATION  OF  MASS  NORMALIZED  MOMENTUM  PER  UNIT  LENGTH 
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Vi,Vj+l  are  the  sPeeds  (i*6*.  magnitude  of 
velocities)  of  projectile  size  j 
and  j+1  respectively 

Xj*Xj+l  are  the  final  displacements  of  the 
debris  particles  from  their  initial 
position. 

If  M  is  multiplied  by  the  mass  of  one  panel  then  the  actual  mo¬ 
mentum  per  unit  length  may  be  obtained.  However,  M  is  presently 
left  in  a  mass  normalized  condition  because  this  allows  for  win¬ 
dow  openings  in  the  panel  and  variation  in  material  properties. 
Presently  H  is  utilized  to  form  three  different  relationships; 

•  Mass  normalized  cumulative  debris  momentum 
per  unit  length  along  the  profile. 

•  Mass  normalized  minimum  debris  momentum  per 
unit  length  along  the  profile. 

•  Mass  normalized  maximum  debris  momentum  per 
unit  length  along  the  profile. 

Once  the  value  of  M  has  been  determined,  it  is  applied 
along  the  length  of  profile  determined  by  Xj  and  x.+^.  In  the 
case  of  cumulative  momentum,  all  the  individual  H  for  all  size 
ranges  and  for  all  floor  heights  are  accumulated  along  x.  This 
gives  an  indication  of  the  amount  of  projectiles  and  their  total 
effect  along  the  debris  profile.  The  minimum  momentum  relation¬ 
ship  along  the  profile  consists  of  the  minimum  momentum  of  any 
single  size  range  at  eacl  x  location  (i.e.,  every  foot)  along  the 
debris  profile.  The  maximum  momentum  is  likewise  the  maximum 
effect  of  any  single  size  range  acting  along  the  debris  profile. 
The  maximum  and  minimum  momentum  relationships  along  the  debris 
profile  establish  hounds  on  the  individual  projectile’s  momentum. 
Whereas  the  maximum  and  minimum  momentum  bounds  give  the  effect 
of  individual  projectiles, the  cumulative  momentum  is  some  indica¬ 
tion  of  the  effect  of  many  projectiles  landing  at  any  one  spot 
along  the  debris  profile.  These  relationships,  when  coupled  with 
available  biological  data  as  to  impact,  are  sufficient  to  esti¬ 
mate  the  casualties  caused  by  flying  structural  debris. 
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4.4  SAMPLE  PROBLEMS 


Two  example  problems  were  run  on  SINBAD  to  Illustrate  the 
system  versatility.  The  results  of  these  problems  are  only  pre¬ 
sented  to  demonstrate  the  problem  solving  capability  of  the  sys- 
t.  -i.  They  are  not  meant  to  illustrate  actual  debris  situations. 

The  first  problem  is  illustrated  in  Fig.  42  and  includes  four 
free-standing  frangible  walls  all  of  the  same  length.  No  shield¬ 
ing  of  one  wall  by  another  is  assumed  to  take  place  since  the 
example  is  designed  to  show  the  superposition  procedure  alone. 

The  time  to  fragmentation  of  all  walls  is  assumed  to  be  zero  as 
is  the  initial  velocity  of  all  fragments.  In  this  problem  all 
calculated  parameters  (i.e.,  debris  profile,  size  distributions 
at  selected  points  in  the  profile,  and  momentum)  were  printed 
first  for  one  wall,  then  two,  three  and  finally  all  four  walls. 

The  input  to  the  problem  is  also  printed  and  both  it  and  the  out¬ 
put  are  displayed  as  Appendix  C.  Plots  of  the  profiles  resulting 
from  the  different  wall  combinations  are  illustrated  in  Fig.  43 
through  46.  All  profile  distances  are  relative  to  the  first  wall 
and  the  remaining  walls  are  located  down  wind  of  the  first  wall. 

It  may  be  seen  that  this  example  illustrates  that  multiple  wall 
configurations  may  be  studied  and  that  the  walls  can  have  ^i^ffer- 
ent  structural  configurations. 

The  second  sample  problem  which  is  independent  of  the 
first  example  demonstrates  how  a  variation  of  parameters  study 
on  the  aerodynamic  properties  of  a  single  brick  may  be  carried 
out  conveniently  with  the  system.  A  free-standing  wall,  40  floors 
at  400  ft  high,  consisting  of  only  a  single  size  particle  (i.e., 
a  masonry  brick  with  nominal  dimensions  of  2-1/4  x  3-3/4  x  8  in.) 
is  exposed  to  a  1  MT  weapon.  The  brick  has  essentially  three 
orientations:  side-on,  face-on  and  end-on.  The  aerodynamic  prop¬ 
erties  of  these  three  orientations  have  been  documented  previously 
(Ref.  14).  Five  separate  cases  were  run  on  SINBAD.  These  included: 

•  A  volume  equivalent  sphere. 

t  Side-on  orientation. 

•  Face-on  orientation. 

•  End-on  orientation. 

•  The  numerical  average  of  cases  2,  3  and  4. 
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The  results  of  the  analysis  are  summarized  in  Fig.  47  through  51 
for  debris  profiles,  Fig.  52  through  56  for  normalized  cumulative 
debris  momentum,  and  Fig.  57  through  61  for  maximum  and  minimum 
bounds  on  normalized  debris. 

It  is  perhaps  interesting  to  note  the  almost  exact  corre¬ 
spondence  between  cases  1  and  5  of  this  problem.  This  is  to  be 
expected  since  an  object  with  the  dimensions  of  a  brick  is  not 
very  different  in  shape  from  a  spherical  object  when  an  average 
orientation  is  assumed.  Larger  objects  with  more  extreme  dimen¬ 
sional  variation  will  probably  not  display  this  similarity. 

The  two  examples  presented  were  to  show  the  versatility  of 
SINBAD.  Thus  it  is  difficult  to  draw  specific  conclusions  as  to 
debris  dispersal  from  these  two  problems.  The  second  example  how¬ 
ever,  does  illustrate  that  maximum  cumulative  debris  occurs  at  the 
same  point  down  range  as  maximum  debris  depth.  This  fact  is  sub¬ 
stantiated  by  Table  7.  Intuitive  reasoning  would  lead  to  this 
same  conclusion  since  the  point  of  maximum  debris  height  is  more 
than  likely  the  point  where  the  most  individual  particles  fall. 

The  maximum  momentum  of  an  individual  particle  falls  much  closer 
to  the  wall  than  the  maximum  cumulative  momentum. 

Table  7 


SUMMARY 

OF  RESULTS 

OF  EXAMPLE  PROBLEM  2 

Particle  Type 

Maximum 

Height 

Distance  @ 
Maximum  Height 

Maximum 

Distance 

Sphere  (Fig.  47) 
Side-on  (Fig.  48) 
Face-on  (Fig.  49) 
End-on  (Fig.  50) 
Average  (Fig.  51) 

0.72 

0.60 

0.52 

1.25 

0.73 

210 

250 

290 

no 

190 

380 

480 

560 

220 

370 

Particle  Type 

Momentum 

Cumulative 

Distance 

Sphere  (Fig.  52) 

26 

210 

Side-on  (Fig.  53) 

22 

250 

Face-on  (Fig.  54) 

19 

290 

End-on  (Fig.  55) 

44 

110 

Average  (Fig.  56) 

27 

190 

DEBRIS  HEIGHT  (FT.  1  (XUrt 

IS.  COO  34.000  33.000  40  000  48  000  S6.000  64.000 


CtST0\r.r  rce*  3P!&1N«1  WALL  IF*. !  1X10*1 

onus  PROFILE  or  SUCK  PARTICLE  IN  EMHM  ORIEHTATION 


MASS  NORMALIZED  CUMULATIVE  MOMENTUM  PER  UNIT  LENGTH  tl/SEC) 


-.QOQ  1.000  2.000  3.000  4.000  5.0CD  „  6.000 

DISTANCE  FROM  ORIGINAL  WALL  (FT.)  iX102) 

Fig.  52  CUMULATIVE  MOMENTUM  ALONG  DEBRIS  PROFILE  FOR  EQUIVALENT 
SPHERICAL  PARTICLE 
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NORMALIZED  CUMULATIVE  MOMENTUM  PER  UNIT  LENGTH  (l/SEC) 


MASS  NORMALIZED  CUMULATIVE  MOMENTUM  PER  UNIT  LENGTH  (l/SEC) 


C3 

O 


Fig.  55  CUMULATIVE  MOMENTUM  ALONG  DEBRIS  PROFILE  FOR  END-ON 
ORIENTATION 
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■  l  P  I  I  F  .  F  I  I 

-.QOQ  1.300  ?.30Q  3.000  4.300  5.C0C  ,  6.000 

DISTPNCF  FROM  ORIGINAL  WRLL  (FT.)  IX! 02) 

Fig.  58  MAXIMUM  AND  MINIMUM  MOMENTUM  ALONG  DEBRIS  PROFILE  FOR 
SIDE-ON  ORIENTATION 
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MASS  NORMALIZED  MAXIMUM  AND  MINIMUM  MOMENTUM  PER  UNIT  LENGTH  CI/SEC) 


END-ON  ORIENTATION 


MASS  NORMALIZED  KnXIMUM  RNO  MINIMUM  MOMENTUM  PER  UNIT  LENGTH 


4.5  FRAGMENTATION  DELAY  TIMES  AND  INITIAL  VELOCITIES 


As  previously  mentioned,  the  trajectory  analysis  used  in 
the  SINBAD  code  is  based  on  the  assumptions  of  zero  initial  vel¬ 
ocity  and  zero  fragmentation  time.  The  trajectory  analysis  is 
essentially  a  numerical  solution  of  a  complex  differential  equa¬ 
tion.  Since  this  numerical  solution  can  have  arbitrary  initial 
conditions  (i.e.,  delay  time  and  initial  velocity),  a  study  was 
made  to  see  how  changes  in  the  fragmentation  delay  times  might 
s! lect  the  final  transport  position  of  a  particle.  The  results 
of  that  study  are  summarized  in  Fig.  62.  The  figure  illustrates 
the  influence  of  delay  time  on  the  final  distance  a  projectile 
travels.  Case  A  is  for  a  particle  initially  at  271  ft  above 
ground  surface,  while  Case  B  is  for  a  particle  at  31  ft  above 
ground  surface.  As  delay  time  is  increased,  the  total  distance 
a  particle  travels  decreases.  This  decrease  however,  is  insig¬ 
nificant  for  delay  times  which  are  physically  meaningful (i.e. , 
up  to  0.1  sec)  for  frangible  panels  commonly  found  in  struc¬ 
tures.  The  delay  time  variation  was  made  again  with  a  zero  in¬ 
itial  velocity.  Increasing  initial  velocity  wilJ  tend  to  off¬ 
set  the  delay  time  effect.  As  fragmentation  time  for  an  element 
increases,  the  strain  energy  within  the  element  builds  up.  This 
strain  energy  is  likely  to  impart  some  kinetic  energy  to  the 
particle  when  it  is  free  to  fly.  Therefore,  an  increase  in  frag¬ 
mentation  delay  lime  tends  to  be  counteracted  by  a  corresponding 
increase  In  initial  velocity  and  t’*e  entire  effect  on  total  par¬ 
ticle  displacement  is  negligible. 

4.o  MODIFICATION*  OF  BLAST  LOADING  Dl*K  TO  IX)CAL  SHIELDING 

One  companion  problem  associated  with  debris  estimation 
is  an  accurate  desript foo  of  the  Mast  loading.  Most  estimates 
of  blast  loading  on  structures  are  developed  under  the  assumption 
that  there  are  no  obstruct  ions  between  ground  zero  and  the  point 
of  load  application.  In  the  real  world  problem  this  is  far  from 
true;  the  blast  wave  must  interact  with  a  variety  of  obstructions 
in  its  path  to  the  structure  of  Interest.  This  phenomenon  is 

III  «*  S  t  A  «C  H  INSfltUlt 


ION 


Fig.  62  INFLUENCE  OF  FRAGMENTATION  DELAY  TIME  ON  FINAL  TRANSPORT  DISTANCE 


known  as  blast  shielding.  Blast  shielding  is  accounted  for  in 
the  SINBAD  model  by  attenuating  the  free- field  overpressure  by  a 
factor  which  is  an  empirical  function  of  building  height,  length, 
and  spacing  from  contiguous  structures.  This  factor  was  deter¬ 
mined  in  a  previous  experimental  program  (Ref.  15)  and  Fig.  63 
illustrates  the  applied  results.  The  three  curves  represent 
different  ratios  of  exposed  length  to  height  for  the  structures 
investigated.  The  separation  ratio  is  determined  from  the  spa¬ 
cing  between  neighboring  structures  and  the  height  of  the  struc¬ 
ture.  These  curves,  Fig.  63,  are  based  on  previous  model  studies 
and  are  the  most  appropriate  data  which  could  be  found  on  atten¬ 
uation  due  to  structural  shielding. 

4.7  IMPINGEMENT  OF  DEBRIS  FROM  ONE  STRUCTURE  ON  ANOTHER 

Although  it  is  conceivable  that  under  the  right  set  of 
circumstances  the  debris  from  one  structure  might  collide  in 
midflight  with  another  structure,  this  result  has  not  been  ob¬ 
served  in  problems  run  to  date.  Such  a  result,  in  any  case,  is 
difficult  to  observe  and  still  more  difficult  to  analyze.  This 
phenomenon  has  not  been  incorporated  into  SINBAD  and  can  only 
be  detected  from  intermediate  results. 

4.8  INTERIOR  BUILDING  CONTENTS  AS  POTENTIAL  DEBRIS 

After  the  blast  wave  interacts  with  the  exterior  walls 
of  a  structure,  it  enters  the  interior  of  the  building.  During 
the  transition  from  the  outside  to  the  inside  of  the  structure 
the  blast  overpressure  undergoes  still  another  attenuate  '  .  This 
attenuation  is,  again,  determined  from  empirical  results  obtained 
from  an  experimental  investigation  (Ref.  16).  As  discussed  pre¬ 
viously,  the  SINBAD  Code  operates  on  an  idealized  space  consis¬ 
ting  of  lumped  particles  at  discrete  initial  heights. 
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These  particles  may  be  any  size;  however,  vt  some  finite  size  the 
model  may  not  yield  plausible  results.  For  interior  items  such  as 
furniture,  ice  boxes,  etc.,  SINBAD  will  yield  good  results,  how¬ 
ever,  large  bulky  objects  such  as  might  be  found  in  a  warehouse 
are  another  story.  These  objects  are  highly  sensitive  to  dif¬ 
fraction  loading  and  must  gain  some  inertia  before-  they  can  be 
picked  up  by  drag  loading.  A  more  meaningful  analysis  for  this 
type  of  interior  debris  item  might  include  a  sliding  overturning 
study  that  would  establish  whether  the  debris  can  start  moving 
or  not.  When  it  is  established  that  the  debris  moves,  the  SINBAD 
analysis  may  be  utilized. 
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C  LIMIT  ANALYSIS  OF  framfs 

rt MANSION  A (20,30) ,S(60>  *ASAT (20,20)  , INDEX <20) ,P(2P> 

rr mens i on  satxoo)  ,pmooj  ,alf oo>  ,cx (2o> 

DIMENSION  CM(30)  ,D<60)  ,H(30) 

c 

c 

C  INPUT  data-  np  is  total  degrees  op  PRFEHOM 
c  nf  is  twice  the  numrer  of  members 

C  A  IS  THE  INTERNAL-EXTERNAL  ^ORCE  MATRIX 

C  P  IS  THE  LOAD  VECTOR 

c  PM  IS  THE  PLASTIC  M0MFM  AT  EACH  NODE 

C  SIS  THE  SLOPE-DELECTION  STIFFNESS  matrix  C I N  VECTOR  FORM) 

C 

c 

1  RFAT  {*>,?)  JJ 

2  FORMAT (IIS) 

IF  (JJ)  4,4,3 

3  RE AC  (5,101)  NP , NF 

101  FOPMAT  ( 2  I S ) 

READ  (5,102)  ( (A (I , J) , JS1,NF) ,I=1,NP) 

102  format  (bfic.A) 

NFT?SNF«2 

RFAC  (5,102)  (S(I) ,I=1,NFT2) 

NH 1 *NP  *  1 

READ  (5,10?)  (P(I),I=1.NP) 

PEAT  (5,102)  (PM( I) , 1=1, NF) 

WRTTF  (6,103) 

103  FORMAT  (31M1LIMIT  ANALYSIS  OF  RIGID  FRAMES//) 

WR I TF (6,5)  J J 

5  F0RmAT(6H1CmSL  NO, I  3) 

WRITE  ( 6 , 1 0 A ) 

10 A  POOMAT  (13M0TMF  MATRIX  A) 

TO  I0«j  U1,NP 

105  WRITE  ( b , 1 06 )  !,(A(I,J),J=1,NF) 

106  FORMAT  (AH  ROW , I  3 ♦ 1 X , IP 4F1 6 , 7/ ( 8X , IP At  1 6 . 7 ) ) 

WRITE  (6,10?) 

107  FORMAT  (13H0THE  MATRIX  S) 

DO  10«  1=1, NF 

11=11-1) /2°Z  +  1 
12=  C I ♦ 1 ) /2«2 
13s2«I-1 
I<«=2»I 

106  WRITE  (6,109)  I,n»S(l3),I?,S(I<0 

109  FORMAT  (AH  ROWi 1 3»5x»3hCOL , I  3, lpF16. 7,Sx,3MC0L, 13, 1PE16, 7) 

WRITE  (6,110) 

110  FORMAT  (13H0THE  MATRIX  P) 

DO  221  I s 1 , NP 

221  WRITE  (6,106)  I,P<!) 

WR I TF  (6,111) 
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r>  r>  nr>n  orto 


111  FO«*ATU*HQTH£  MATRIX  PM) 
!30  171  1*1. NP 
7ft  WRjTt  (6.10*)  I.PmU) 

CO  112  Isl.NF 

c 

c 


INVERT  THE  STIFFNESS  MATR 1 X  .  D«<$**«1> 


TF  <I/2«2-I>  113.112.112 
113  Nl*2*l-1 
N?s2«I 
IP1=I*1 
N3s?*IPl-l 
N«s2*IPl 

IF  ( S ( M ) )  711.712.711 

712  D(N«)8l./5(N4> 

n(Nl)=r>. 

0 ( N2) aO • 

D (N3) a0» 

00  TO  112 

711  IF(S(N4)>  710.713.710 

713  D(Ni>al./S(Nl) 

P(N2)sO. 

H (N3) aO, 

r(N4)sO. 

BO  To  112 

710  TFmPsi./(S(n1)»S(N<*)-S(N2)#S(N3)  ) 
T ( N 1 ) sS ( N4) «T£MP 
C(N4)sS(Nl)«TEMP 
C(N2)a-S(N2)  OEMP 
f ( S3) *0 ( N2 ) 

11?  CONTINUE 


FL AST  1C  ANALYSIS  FOR  UNIT  LOADS 


500  NCYCLSO 
Cl F*0, 

DO  *4  1*1 »NP 
24  CxtDsO, 

rn  U  i*i. nf 

26  Cm  1 1 >  so • 

15  ro  u*  I *1 . P  P 
m  116  Jal »Np 

ASAT(l,J)sn, 

TO  lib  x*l.Nr 

K?»<*»1 )/?•? 

X  3«?<*K-1 
Xbs2«« 
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UUUULiU  C-UU  U 


116  ASAT  ( i  ,  j)  =ASAT  ( I  *  J)  +A  ( I  ,K  )«  C5  C K  3 >  «A  ( J,K1)  +5>  (K4)  <*A  ( J»K2)  ) 
CO  151  1=1, NP 

151  ASAT  ( I , NP 1 ) *P { I ) 

DC  U?  1=1, NP 

117  !«rfX(I)sO 
116  AmAXs-i. 

CO  119  I=1,NP 
T  F  <INI?CX<I)>  119,120,119 
120  TtMPsAP.S(ASAT(!,I)) 

T  F  (TENP-AMAX)  119*119,121 


121  IR0W=I 

AKA*=TEKP 
119  CONTINUE 


CHECK  FOR  ZERO  IN  PIVOT  ELEMENT  OR  EXCESSIVE  DEFLECTION 
TnF.SE  CONDITIONS  APE  INSTABILITY  CHECKS 


IF  UMAX)  122,147,124 

124  I V  T  F  X  (IPOW)  =1 

PI VOTsl ./ (ASAT (IPOW, I  ROW)  ) 

DO  125  J s 1 ,  N P 1 

125  ASAT ( IPOW, J) =ASAT (IROW, J)«PIVOT 
n0  126  1=1, NP 

IF(I-IROW)  127,126,127 
1 ?7  TFMp  =  ASAT ( I , I  ROW) 
rc*  126  J  =  1  , N P  1 

128  ASAT ( I , J) sASAT ( J , J J - A S A T ( I POW , J ) «TEMp 

126  CCNT 1 MJE 

GO  TO  !  1 8 

147  WRITE  <*.347> 

347  POPMA7  (24H0ZFPO  PJVOT  IN  INVERSION) 

GO  TO  47 

I??  no  3i i  :=i*np 

IF  UR$lASATU,hPl>  )-l.E*10)  311,647,647 
311  CONTINUE 
GO  to  303 
647  w tt  T  r  F  (6,647) 

647  FORMAT  ( 21H0DEFL  EC  T ION  TOO  LARGE) 

GO  T0  47 


COMPUTE  THE  von ENTS 

303  DO  T  ? 1  1  =  1, NF 
Il=< I-l ) /2«2*1 
!?=<!♦! ) /?*2 
I3=*"I-1 

!4=?»I 
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u uu uu  ou 


5ATX(!)*0. 

PO  131  *xl,NP 

131  $1TX(  I)«SAT*  (n^ASAT  IK  n  |«A(K  ,!!)♦$(  I*)«AfK«|2|) 


F I ND  ADDITIONAL  LOAf  F AC  TOR  PE3UI0FD  TO  MING  JOINT  WITH  LARGEST 
mOM£NT  UNt1E.fi  UNIT  LOAD  UP  TO  PlASTjC  M0MENT 

check  that  moment  is  increasing  **  *agnitupe  under  unit  load 


mo  201  I*1*NE 

TF  (ARS(SATKd)  )-l.E-0*)  202,202,203 
2n2  AlE<II*1.E?0 
GO  To  201 

203  ALF<l)»<PM<n-ARS<CM(l)  n/A6S(SATX(I>) 
201  CONTINUE 


3AlFei,E20 
m  204  1*1 *NF 
test  s  C*MI)»SATX<I) 

IF  (TEST)  ?34,?0S,20?> 

205  IFCALEI  n-SALf  )  1206.204,204 
1206  SAIF  a  ALF(I) 

NEH=I 

204  continue 
c 
c 

c  if  there  is  no  increase  in  loat  factor, collapse  mechanism  exists 

c 

c 

IF  ISALE-l.f-OT)  247,247,302 
247  WRITE  16,447) 

447  FORMAT (2?HPlOAr  FACTOR  TOO  SMALL) 

GO  7C  4? 

C 

C 

c  commute  moments  undfr  current  load  factor 

c 

c 

302  Pr  POT  |«l »\F 

SAT«U)*SALF*SATX(I) 

207  Cm(I)«CmU)*SATx(H 

r 

c 

C  COUMLE-CmECk  An**lSSiOlMTT  ex  SOluTrON.“CM|NT,LE.*M  at  ALL  jC (NTS 

c 

c 

Po  *1*  I*1*nf 

if  fPMtn«ApS<CMt!n»],E*03)  547,314,314 
314  CONTINUE 
*f  304 

547  *s  f  Tr  14,7471 

747  FOR*' A T (24!*m»l AST  1C  MOMfNT  fxCtEOCP) 

On  to  <,* 


C 

c 
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I 

I 

I 

I 

I 

I 

I 


c  WRITE  MOMENTS  AMO  DEFLECTIONS  UNDER  CURRENT  LOAO  FACTOR 

C 

C 

304  ClFsClF+5AIF 
DO  206  I=1»NP 

A  S  A  T  (I«NPi)sS  ALF°A.S  AT  <  X  «  NP  J  > 

20 ft  CX ( t ) =C X ( I >  * A$A T ( I »NP  1 ) 

NCVCl=NCTCL*l 

WRITE  (6.401>  NCYCL.NPw 

401  FORMAT  (16H1PLASTIC  wjnGF  VO . , J 3 , 2 X » 1 5HF0RmED  AT  POINT, 13) 
wRirr  (6,402) 

402  FoRMATU?HCiOAP  FACTOR, 3X,10HAODITIonAL«9X,10HCUMULATIVEJ 
write  (6,403)  NCYCL.SAlF ,CLr 

403  FORMAT ( TWO  STAGE  <,I3»lH),!PEia.7,l»E19,7> 

WRITE  (a, 404) 

404  FORMAT  an4ODEFLfCTIOV,4X,10NADDITIONAL,9X,10MClJMULATIVE/) 
CO  208  1=1 »Np 


208  WRITE  ( 6 , 4 OS ?  NASATU.NRi)  ,CXU) 

405  F0RMAT(3W  X  (,I3,1H)  ,IPE22,7,1PE19.T) 

WRITE  (6,406) 

406  FORMAT  (  ?HLiw0MENT8X,i0HAnnirl0NAl  ,9X10NCUMULATIVE10X,8HPLAS  MOM/) 
DO  209  I  =  1  * NF 

209  WRITE  (6,407)  I , S A T X ( I ) ,CMfI),PM(l) 

407  FoRMATOH  M(,I3,lw>  ,F  lfl.4,?F19.4) 

C 

C 

c  COMPUTE  INELASTIC  hinge  ROTATIONS  UNDER  CURRENT  LOAD  FACTOR 

C 

c 

DO  933  1  =  1 , N E 
11  =  0-1) /J!*2*l 
12=  U ♦  1 )  /?.»Z 
I3=2«I-1 
T4=2«I 

Hr?)sD(I3)PCM(Il)+P(I4)«CM(I2) 

DO  933  Xsl.NP 

933  H(I>=H(I)-A(X,I)OCX(X) 

DO  501  1=1, NF 

501  IF (APS (H( I) ) ,LT. (l.f-07) )H( I)s0,0 
WRITE  (6,13*) 

138  Format (1H0,14X,15HHINGE  ROTATIONS/) 
m  939  1=1, NF 
939  WRJTF  (6,140)  I,H(I) 

140  format ( j oh  at  point r,i3, im) ,ipe15*7> 

c 

c 
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MODlfy  STIFFNESS  MATRIX  TO  INCLUDE  PLASTIC  HINGE 


IF  <fNPH/2«2>-NPH)  211,210,210 
211  N1s2®NPH-1 
N'2s?«NPH 
NPHPlsNPH+l 
M3s2«NPhPl-i 
N4s2«NPnPl 

S(NH)sS(N^)«a.-S(N2)«S<M3)/(S(»l  1I«S(N4)  )) 

S  ( V  1 )  s  0  • 

S (N?)=0« 

S (N3) sO, 

Gn  TO  ?12 
210  nphw1snph-1 
m=2«nphmi-i 

N'2s2®NPMM1 

N3s2#NPh-1 

N4s2®NPH 

S(NI)s$(N1)#(1.-S(N2)«s(N3)/{S(N1)«S(N4M) 

SIN?) sc# 

S  IN T ) sO. 

S  IN4) so, 

212  GO  TO  lb 
C 

c 

47  WRJTF 


FORMAT(36HOCOLLAPSE  MECHANISM 

GO  TO  1 
4  STO° 

ENT 


HAS  BEEN  RfACHFD) 
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APPENDIX  B 


RESULTS  OF  SAMPLE  PROBLEM 


CASE  NO  2 


THE 

ROW 


MATRIX  A 


ROW  2 


ROW  3 


ROW  4 


ROW  5 


ROW  6 


ROW  7 


ROW  ft 


-0  .GC-nOOOOE 
-P.  OijO'.'COOE 

•  0 .  0  Q  0  Ci  0  0  0  E 
1 . 0'J'VPC'OOE 
-O.OO.u.OGOE 
-O.DCKU'OOCE 

-  P  .  0  0  r'  i’  0  0  (i  E 

-C  .  C'j 'cOPOE 
-O.nQ.'ftOrjCE 

-  G .  C  0  •}  >)  n  0  0  £ 
-C.OO.V.'OOOE 
-O.0i;-)(.>A0OE 
-O.OOOuOnoE 

-C .  iIpnoOuGE 

-C .  nnor.nnoE 

-n . ooooouoe 
"O.OOOC-OCOS 
-0.0  fi  0 £  j  o  E 
-0 . 0 1 )  i 0  ».i  0  E 
l.OjO’jO'KE 
-0 t  OoncOCGE 

"0 . 0  0  f:  .,*?•:  OE 
- G .  0  J 0 tiP o 0  E 


-39 
-39 
-39 
00 
-39 
-39 
*39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
-39 
00 
-39 
-39 
-39 

-fi,ocr.i;.r»or.E-39 
-O.OonrOuOE-39 
-0,OCijt.)OOOE-39 
-  0 . 0  0  0  0  0  0  0  E  -  3  9 
-f',OnooonOE-39 
-0 , 0(J'’-:(i,;>oe -39 

-0 . 0  fi  I-: , ,  n  <  >  C  E  -  3  9 
1 .0.J'JCOOOE-Ol 
I  *  0 1) : 1 C’  Ci  OOE-Ol 

1 ,9yr.rnoCE-oi 

JOE-39 

-o  .0ooC'0f.'0E-39 

•v»  00,. [.CODE -39 
-0 ,0ijfi*.C00E-39 
-O»O0':cO0Oe-39 


l.OOOOOOOE  00 
•  nnO000CE-39 
-0 • OOGOOOQE-39 

-0 .  oonooooE-39 

-  0  •  00  0  0  0  0  0  E  -,3  9_ 
-0  •  0  0  i>  0  0  0  0  E  -  3  9 

i.noOOOOOE  00 

-  0 .  <1 0  0  0  0  0  0  E  -  3  9 
l.nuOOOuOE  00 

-C.oooononE-39 

-O.OOOpOOOE-39 
-o.onOuoooE-39 
l.ntjoooonE  00 

-  0 . 0  0  0  0  0  0  0  E  -  3  9 
-0. nOOOOOOE-39 

•  0 . 0  0  0  0  0  0  0  E  -  3  9 
-o.  iinoooorjE-39 

-0 .  D,j0yO0r'E-39 
-O.OOJOOOOE-39 

•  0  •  00-',0900E-39 
-0. nOOOOO'JE-39 
-0 .  OOOOOO'iE-39 
-O.noOOOOoe-39 
-0. nOOOOOOE-39 

l.OUOOOOOE  00 
"0  .OOOOOOOE-39 
■  0  •  (1  j  0  0  0  0  0  E  -  3  9 
•0.  nOOOOOOE-39 
•n •  OUOuOOuE-39 
.0.!i00000oe-39 
1  ..-ionuOOOE-01 
1*  nnoooooE-oi 
l  .nnoooooE-oi 

-0  .tiOOOOOOC-39 
-0  .r,000P00E-39 
■0  .nnC00u0E-39 
O0Oon0rE«39 
•0  .'iononOCE-39 
■U, OOOOOOOE-39 


l.OOOOOOOE  00  -0. 

-0. OOOOOOOE-39  -0. 
-0. OOOOOOftE-39  -0. 
-0. 0000O00E-39  -0. 

-O.OOOOOOOE-39_  -0. 
-0. OOOOOOOE-39  -0, 
l.OOOOOOOE  00  -0. 

-0. OOOOOOOE-39  -0, 
l.OOOOOOOE  00  -0. 

-O.OQOOOnOE-39  -0. 
-n. oooooooe-39  -Oj 
-  0 • 00  0  0  ^  0  0  E - 39  -0. 

l.OOOOOOOE  00  -o. 

-0. OOOOOOOE-39  1. 
-0. OOOOOOOE-39  -0. 
-O.OtiOOOOOE-39  l. 
-o.  oonoonoE-39  -o. 
-0.0n00O00E-39  -o. 

-0. OOOOOOOE-39  — 0 • 
-0.  OOOOOOOE-39  -o. 
-0. OOOOOOOE-39  -0. 
-0. OOOOOOOE-39  1* 
-0. OOOOOOOE-39  -0, 
-o.O(jnonbOE-39  -0. 
l.OOOOOOOE  00  -o. 

-0. OOOOOOOE-39  -0. 
-O.OOOOOOOE-39  -0. 
-0. 0000O00E-39  1. 

-0.0000Oooe:-39  -0. 

-O.OOOOOOOE-39  1. 
•i  .ooooonoE-oi  -l. 
•  1  • OOOOOOOE-0 1  -1» 

•1 .000000OE-01  -1. 

-o.onooonoE-39  -o, 
.O.OOOOOOOE-39  -o, 
1  .OOOOOOOE-Ol  3, 

i. ooooonoE-oi  l. 

1  .OOOOOOOE-Ol  1. 
•0. OOOOOOOE-39  -0. 


-  '  .  0  0  o 

i.O 

L.oE-39  -0*  jononocE- 

39  - 

o. 

OOOOnonE-39 

-1» 

the 

ROW 

MATRIX  s 

1  C°L 

1 

7.2oncOunE 

05 

COL 

2 

3. 60000nnE 

05 

! 

ROW 

2 

COL 

I 

3.6r ooCCpE 

05 

COL 

2 

7.2rt00000E 

05 

I 

ROW 

3 

COl 

3 

7.2l'OOOOoE 

05 

COL 

4 

3.6onooncF 

05 

ROW 

4 

COL 

3 

3.  t<  OOOOOE 

05 

COL 

4 

7.200n000E 

05 

! 

ROW 

5 

COL 

5 

7 . 2  C  0  0  0  C  o  e 

05 

COL 

6 

3.6nooonoE 

05 

i 

1 

ROW 

6 

col 

5 

3.6POOOOOE 

05 

COL 

6 

7;2nonooOE 

05 

ROW 

7 

COL 

7 

7 . 2  '*0000oE 

05 

COL 

6 

3-  snooonoE 

05 

] 

ROW 

8 

COL 

7 

a.fei'OOOUOE 

05 

COL 

8 

7,290nonoE 

05 

ononoooE-39 

OOOOOOOE-39 

OOOOOOOE-39  ! 

OOOOOOOE-39  \ 

OOOOOOOE-39 

ononoooE-39  i 

OOOOOOOE-39 

OOOOOOOE-39 

OOOOOOOE-39  i 

OOOOOOOE-39 

OOOoOOOE-39 

OOOOOOOE-39 

0000000t”39 

OnOoOOOE  00 

OOOOOOOE-3V 

OOOOOOOE  00 

OOOOOOOE-39 

OOOOOOOE-39 

OOOOOOOE-39 

nn00000E-39 

OOOOOOOE-39  j 

OOOOOOOE  00 

nnOtioooE-39 

0noon00E-39 

oooooooe-39  j 

OoOOOOOE-39 

OOOOOOOE-39 

ooonoonE  00 
OOOOOOOE-39 
OOOOOOOE  00 
onoooooE-oi 
noonoooE-ol 

OOOOOOOE-Ol 
0.0000  00  E -39 
ooonnot)E-39 
OOOOOOOE-01 
OoOOOOOE-O 1 
OOOOOOOE-Ol 
OOOOOOOE-39 
OOOOOOOE-39 
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ROW 

9 

COL 

9  7 

• 2000000E 

ROW 

10 

COl 

9  3 

^6000U00E 

ROW 

11 

cQl 

11  7 

,2uoooonE 

ROW 

12 

col 

11  3 

Z 60  00  00  0  E 

ROW 

13 

COL 

13  3 

.6000C00E 

ROW 

19 

COL 

13  1 

,6000000E 

ROW 

15 

COL 

15  3 

.60000CHE 

ROW 

16 

COL 

15  1 

,8w0GU0f*E 

ROW 

17 

COL 

17  3 

#  6  •'  0  o  o  c  n  e 

ROW 

18 

COL 

17  1 

.eoouoooE 

ROW 

19 

COl 

19  3 

.6u0000nE 

ROW 

20 

COL 

19  1 

.  8  C  0  u  0  0  ft  E 

THE 

MATRIX  P 

ROW 

1 

-  0 , r  0  ? 

r  o  o  o  E  • 

39 

ROW 

2 

•0,000 

rCuCE- 

39 

ROW 

3 

-u.fle*' 

OOOOE- 

39 

ROW 

9 

-  n .  0  u  i- 

rP'-.CE- 

39 

ROW 

5 

•0 • ooo 

i'-vCOE- 

39 

ROW 

6 

-C.Oon 

i 1 0  u  0  E  - 

39 

ROW 

7 

-".OGI* 

oOCf'E- 

39 

ROW 

6 

1 . n  0  * 

..r-joE 

00 

THE 

MATRIX  PM 

ROW 

1 

6.  n0;i 

voooe 

03 

ROW 

2 

6.0CC 

it’OCOE 

03 

ROW 

3 

6.00' 

.  ,i-?r*E 

03 

ROW 

9 

6. our 

OPOOE 

03 

ROW 

5 

6 .  '+>  u  t 

,  •*'  0  o  0  E 

03 

ROW 

6 

6.00*' 

*  «.•  0  u  0  c 

03 

ROW 

7 

6,  Oor 

.  i .  *  r>  n  £ 

03 

ROW 

8 

6.00- 

<  sOO^t 

03 

05 

COL 

10 

3.6000000E 

05 

05 

COL 

10 

7.2000000E 

05 

05 

COL 

12 

3.6000000E 

05 

63 

COL 

12 

7.2000000E  05 

05 

COL 

19 

i.enoooooE 

05 

05 

COL 

19 

3 , 6000000E 

05 

05 

COL 

16 

i.eoooonoE 

05 

05 

COL 

16 

3.6noonnr>E 

05 

05 

COL 

18 

1.8000000E 

05 

05 

COL 

18 

3»  60000nnE 

05 

05 

COL 

20 

1.8000000E 

05 

05 

COL 

20 

3.60000HOE 

05 

PLASTIC  HINGE  NO*  1  fORMED 
LOAD  FACTOR  ADDITIONAL 
STAGE  I  1)  2.3983869E  03 


AT  POINT  6 
CUMULATIVE 
2'3*63S69£ 


03 


deflection 

additional 

x( 

1) 

-9.2293903E-03 

x< 

2 ) 

-7.1689S83E-03 

x< 

3) 

-9.2293903E-03 

X( 

9) 

-7.2281959E-03 

Xt 

5) 

-9.7192351E-03 

XI 

6) 

-7.22819571-03 

X( 

7} 

7.8952407E-02 

X( 

6) 

1.8986959E-01 

cumulative 

-9.2293903E-03 
-7.1689583E-03 
-9.22939Q3E-03 
-7.228195AE-03 
-4.7192351E-03 
-7.2281957E-03 
7 . 6952607^-02 
1.89869591-01 
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MOMENT 

ADDITIONAL 

CUMULATIVE 

PLAS  MOM 

M< 

1) 

5204.3008 

5204.3008 

60on.oonn 

M< 

2) 

1681.7204 

1861.7204 

6000,0000 

M< 

3) 

2731.1828 

2731.162J 

6000,0000 

M* 

4) 

3451.6130 

3451.6130 

6000. O000 

M( 

5) 

5946.2365 

5946.2365 

6000.0000 

M  ( 

6) 

3365.5915 

3365.5915 

6000.0000 

M< 

7) 

5118.2797 

5118.2797 

6000,0000 

M( 

8) 

5999.9999 

5999.9999 

6000,0000 

M< 

9) 

5204.3008 

5204,3008 

6000. 0000 

M< 

10) 

1881.7204 

1881.7204 

6000.0000 

M  ( 

11) 

2 731.1828 

2 731.1828 

6000.0000 

M  ( 

12) 

3451.6128 

3451 .6128 

6000.0000 

M< 

13) 

-4612.9030 

-4612.9030 

6000.0000 

M< 

14) 

-4241.9352 

-4241.9352 

6000.0000 

M< 

15) 

-4241.9352 

-4241.9352 

6000.0000 

M( 

16) 

-4612.9030 

-4612.9030 

6000. noon 

M< 

17) 

-3451-6126 

-3451.6126 

6000.0000 

M  ( 

18) 

-2999.9998 

-2999.9998 

6000.0000 

M< 

19) 

-2999.9998 

-2999.9998 

6ooo. ooon 

M* 

20) 

-3451.6127 

-3451.6127 

6000. rooo 

hinge  rotations 

AT 

POINT ( 

1» 

0.000u"V,E-39 

at 

POINT ( 

2’ 

0. 0A',<V.i)oE-39 

AT 

POINT ( 

3» 

■j.OCCOr.Df  £-39 

AT 

POINT) 

4) 

•..0Ccr><'ntE-39 

AT 

POINT) 

5> 

OuCOO‘V,E-39 

at 

POINT* 

6' 

.  . uOCO''" JE-39 

at 

POINT) 

7) 

t  .:J^CCfnL C»39 

AT 

POINT) 

8* 

'>•  nruf.r  jiE-39 

at 

point) 

9) 

i  0  G  t  *  0 .7  C  -  3  9 

at 

POINT) 

10) 

“.•iOCorn’jE-39 

AT 

POINT) 

11) 

u0c'i.'P')uE*39 

AT 

POINT* 

12) 

o.ocuOon  E-39 

AT 

POINT) 

13) 

•!.  J0<.H'n'jiE-39 

AT 

POINT \ 

14) 

» ;  •  0  0  0  c  r  of-39 

AT 

POINT* 

15) 

(.•U0'J0n,M.E-39 

AT 

POINT* 

16> 

i’»00uf.r'1uE-39 

AT 

POINT* 

17) 

o.jeuo<‘'^E«39 

AT 

POINT) 

18) 

.;.:>rtC.'’PLF-39 

AT 

POINT) 

19) 

f»0 (»'’,! CE-39 

AT 

POINT) 

20) 

•  I.CCl'n,':’,  F-J9 

PLASTIC  HINGE 

no*  2  formed  at 

POINT  5 

LOAD  FACTOR 

additional 

cumulative 

STAGE  I  2) 

1.8046505E  01 

2.36643341  03 

deflection 

additional 

cumulative 

XI 

1) 

-9t4618827£-05 

-9.3240091C-03 

XI 

2) 

-5.76485S4E-05 

-7 .22610685-03 

XI 

3) 

-9*46186292-05 

-9.324009)1-03 
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x< 

4) 

#■ 

•  I «  1446163E*  04 

X* 

51 

5.7230812E~^ 

x< 

6) 

m 

>1 , 1446162E-04 

xt 

7) 

6«899725U-04  . 

X( 

6) 

2.0546C02E-03 

MOMENT 

additional 

m( 

1) 

40.4543 

Mt 

2) 

6.3915 

M  ( 

3) 

38.0481 

Mt 

4) 

30.9046 

Mt 

5) 

53.7635 

M( 

6) 

33.0101 

M< 

7) 

42.5597 

Mt 

8) 

0.0000 

M< 

9> 

40.4543 

Mt 

10) 

6.3915 

Mt 

11) 

36.0480 

M< 

12) 

30.9046 

M< 

13) 

-44.4395 

M< 

14) 

-37.7849 

Mt 

15) 

-37.7849 

Mt 

16) 

-44.4395 

M< 

17) 

-30.9046 

Mt 

18) 

-0*00  0 1. 

Mt 

19) 

-C.  00.1'. 

M< 

20) 

-30.9046 

hinge  rotations 

AT 

POINT  < 

1) 

.  niiOG E-39 

AT 

POINT ( 

2) 

(;*0OO0r9«.  E-39 

AT 

POINT* 

35 

O.0000n«n..E-39 

AT 

POINT  < 

4) 

u.i'COC'*’)  -E-39 

AT 

POINT  f 

5» 

*. .  fjOuorof.  E-39 

AT 

POINT ( 

6> 

1  .uOOoic'tE-39 

AT 

POINTt 

T) 

i,.  joor**.."‘g-39 

AT 

POINTt 

6> 

-2.3310143C-04 

AT 

POINTt 

91 

f:  •  GOOO^uoE-39 

AT 

POINT* 

10) 

« r  •  *  *  f)  G  r  ■  j  E  —  3  9 

AT 

POINTt 

ID 

.  ;}Cvr  ^  .1  ,E-39 

AT 

POINTt 

121 

,  *uOi  j-'  ">■  E-39 

AT 

POINTt 

13) 

'•,00'..'  ^ 'E-39 

at 

POINTt 

14) 

-E-39 

AT 

POINT* 

151 

i',.n,0.’>p'E-39 

AT 

POINT* 

16) 

ft.-ji)'.’**.-  •*.  E-39 

AT 

POINT* 

17) 

•  •  ,f:7  >r>  v, .  E-39 

AT 

POINT  t 

10) 

•  .*.*•'  i»o r'  -* 1  ■  E  *  39 

AT 

POINT* 

19) 

t;.ynOC',f)'|-39 

AT 

POINT* 

20) 

..••OCOOOjC-l-St 

-7.3426570E-03 
-4.6620G43E-03 
-7 • 3426573E-Q3 

7,9642S7H-02  . . 

1.9191919E-01 

cumulative  *»las  mom 

5244,7551  6Gno,noon 

1888.1118  6000.0000 

2769.2308  6r0r .0000 

3482.5177  6000.0000 

5999.9999  6000, 0000 

3398.6015  6000. Oooo 

5160,8394  6000.0000 

5999.9999  6000,0000 

5244,755.1  6000.0000 

1888.1118  6000.0000 

2769.2308  6000.0000 

3482,5174  6000.0000 

-4657.3425  6n00.000ft 

-4279.7200  6000,0000 

-4279.7200  6000,0000 

-4657.3425  6000,0000 

-3482.5172  6000,0000 

-2999.9998  6000,0000 

-2999.9998  6000. oooo 

-3482.5174  6000, 0000 


PLASTIC  HINGE  HO*  3  FORMED  AT  POI^r  9 


LOAD  FACTOR  ADDITIONAL . CUMULATIVE 

STAGE  1  3)  2.2768410E  02  2.5941175E  03 


DEFLECTION  ADDITIONAL 
X (  l)  -1 .46O3O46E-03 

X (  2)  -6. 1703011 E“ 04 

X (  3)  -1 .4603046E-03 

X (  4)  -1 .4808722E-03 

X (  5)  7.4043609E-04 

X (  6)  -1.4808722E-03 

X (  7)  1.1860691E-02 

X (  8>  3.0303036E-02 


CUMULATIVE 
-1.07843141-02 
-7.8431369E-03 
-1.0784314E-02 
-8 « 8235291 F-03 
-3.9215682E-C3 
-8 • 8235295E-03 
9  •  1503260E-02 
2 . 2222223E-01 


MOMENT 

additional 

M< 

1) 

755.2449 

M  * 

2) 

229.5353 

M( 

3) 

407.2399 

M* 

4) 

399.8356 

M  ( 

5) 

-0.C*'0( 

M< 

6) 

307,2010 

M  ( 

7) 

662.6904 

Ml 

6) 

Ml 

9) 

755-2449 

Ml 

10) 

229.5353 

M( 

ID 

407.2400 

Ml 

12) 

399.8356 

Ml 

13) 

-636.775 1 

Ml 

14) 

-484.9857 

Ml 

15) 

-484.9057 

Ml 

16) 

-636*7751 

Ml 

17) 

-399-8355 

Ml 

18) 

— 0 . 0  0  0  0 

Ml 

19) 

-  0  •  G  0  0  0 

Ml 

20) 

-399.8355 

CUMULATIVE 

5999.9999 
2117.6471 
3176.4707 
3882.353 2 
5999"  •  9999 
3705.8825 
5823.5297 

5999.9999 

5999.9999 
2117.6471 
3176.4707 
3882.3530 
-5294.1176 
-4764.7057 
-4764.7057 
-5294.1175 
-3882.3527 
-2999.9998 
-2999.9998 
-3882,3528 


AT  POINT (  1> 
AT  POINT  <  2> 
AT  POINT  (  3) 
AT  POINT(  41 
AT  POlNTt  5> 
AT  POINT (  6J 
AT  POINT*  7) 

AT  POlNT(  01 
AT  POINT (  9) 
AT  POINT*  101 
AT  POINT!  11) 
AT  POINT*  12) 
AT  POINT*  13) 
AT  POINT*  14) 
AT  POINT (  15) 
AT  POINT*  16) 
AT  POINT*  17) 
AT  POINT*  18> 
AT  POINT*  19) 
AT  POINT*  20» 


hinge  ROTATIONS 

...  !>(  ,  CfM  .E-39 
.  d  0 1. 0  0  o  ii  E  -  3  9 
oCCHrOi  E-39 
.  0  0  !.■  1?  f'  i..  1 1  E  -  3  9 

.1 .4705887E-03 

<>  uOuCrr';  E-39 
r  .  0  0 :  i  T/  E  -  3  9 

.3.431 3742E-03 

.1.  E-39 

1.  .UOMirG  jf-39 
n  .COl  0 mi: i.. e-39 

•  i .  or C0rf'CE-39 
0  •  uOooncoE-39 
M.000Cr'0-:.E-39 
■ ;  .  ii  :> O'  ;.<E-  39 
{, .  0-^!, E-39 

•  .  DCuOii-: .  E-39 

jCOGoCuE-39 
r''roroi:E-39 
>' .  Un0C'M!VE-39 


PLAS  MOM 

eoon.oooo 
6000.0000 
6000.0000 
600" .oono 
6000.0000 
6^00 .0000 
60CC • 0000 
6000.0000 
60G0 .0000 
6000^*5000 
600r  .0000 
6000 .0000 
6ooo. oono 
6000.0000 
6000.0000 
6000^0000 
6000. 0000 
6000.0000 
4000.0000 
6OC0.O0OO 
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'  *  /-*>. 

*\ 

••  '•tH  .1  •til,: 

plastic  Hinge 

Ho.  4  forced  at*  point  1 

i 

load  FACTOR 

ADDITIONAL 

CUMULATIVE 

* 

1 

f 

STAGE {  4) 

1  • 2065j 1 3C-05 

2.594U75E  03 

DEFLECTION 

ADDITIONAL 

CUMULATIVE 

] 

x< 

1) 

-9.5428928E-H 

-1.0784314E-02 

i  2 

X( 

2) 

-4.0570425E-H 

-7 , 843 1369E-03 

x< 

3) 

-7* 5083876^*11 

-1.0764314E-02 

x< 

A) 

-7.8787207E-11 

-6.8235291E-03 

A 

x< 

5) 

4.1089025E-11 

-3.9215681E-03 

X( 

6) 

-S.5568894E-U 

-6.8235295E-03 

i 

X( 

7) 

8.S323677E-10 

9.1503269E-02 

■  \ 

;  | 

x< 

8) 

1.9207236E-09 

2.2222223E-01 

MOMENT 

additional 

CUMULATIVE 

PLAS  MOM 

t  4 

M( 

1) 

f!  .  0  '1  *  1 

5999.9999 

6000,0000 

•  * 

M  ( 

2) 

o.onoo 

2117.6471 

6000,0000 

M  ( 

3) 

0 .  Or  >n. 

3176.4707 

6*0". *000 

i 

M( 

4) 

o.o  **<• 

3862.3532 

600* .0000 

i 

.  * 

M  I 

5) 

-  r.  •  o  -i  <t  •  i 

5999.9999 

600n.00nc 

M< 

6) 

f  .'J-'liU 

3705.6825 

6000.0000 

M( 

7) 

p  .  0  rt  f 

5623.5297 

6*00.0000 

i 

M  ( 

8) 

0 .  Cr-Ciit 

5999.9999 

6*00.0000 

M  ( 

9) 

-0 .  QnOi 

5999,9999 

6*0*. 0000 

.  . 

M( 

10) 

c .  c  *  n 

2117. 6471 

6000. *000 

M< 

11) 

(< .  0  r  t.i  0 

3176.4707 

600*. 0000 

M( 

12) 

n.c^O'f 

3882.3530 

6*00.0000 

M< 

13) 

-  0 .  •} n  '1  u 

-5294.1176 

6*00. OOOC 

M< 

14) 

-0 ,0^00 

-4764.7057 

6000.0000 

Ml 

15) 

-0,0000 

-4764.7057 

6000.0000 

M< 

16) 

*  0  •  c  o  c  i 

-5294.1175 

6*00.0000 

Ml 

17) 

-C- .  cnoo 

-3882.3527 

6000.0000 

Ml 

18) 

0  .  C  o  •* r 

-2999.9997 

6*00,0000 

Ml 

19) 

-f.  -  OcOi.t 

-2999.9998 

600*. *000 

Ml 

20) 

-o.r,',  of.* 

-3882,3528 

6*nr.  oonr. 

HINGE  ROTATIONS 

AT 

point ( 

1) 

.  or.r^nottE-39 

AT 

POINT ( 

2) 

■  ) . *  *  sj  C  n  u  C  E — 3  9 

AT 

POINT ( 

3) 

0.  •)•..) CC*U*E- 39 

AT 

POINT! 

4) 

•:  .OCOO*-G'fE-39 

AT 

POINT ( 

5) 

.1 .4705  887  E-03 

AT 

POINT! 

6 ) 

O.ijCCP**  'JE-39 

AT 

POINT! 

7) 

'>,unC0',r’uE-39 

AT 

POINT ( 

8) 

-3.4313744E-03 

AT 

POINT ( 

9» 

r  .oco(*oogE-39 

.! 

AT 

POINT! 

10) 

G • 0  C  u  0 T  0  OE-39 

• 

at 

POINT! 

ID 

iJ.CGCGOuOE-39 

•  t 

AT 

POINT ! 

12) 

n .  n0u0n0'*E-39 

1 

AT 

POINT! 

13) 

r, ,  orc-Onr.'"'E-39 

AT 

POINT! 

14) 

■,;.OOL*C'»ri!  E-39 

*  * 

AT 

POINT ( 

15> 

f.»  OI'-jIm'g  'E-39 

j 

AT 

POINT' 

16) 

,,0**uE-39 

•A 

[ 

AT 

POINT ( 

17) 

*.  o  «  o  c  g  -  39 . 

i 

AT 

POINT! 

18) 

i:  •  wocc*  6oE”39 

T 

AT 

POINT  ! 

19) 

o.0*0QOCCE-39 

l 

"A 

AT 

POINT! 

20) 

r  .ooccno  IE-39 

f 
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1 

:4 

1 

* 

PLASTIC  HINGE  NO.  5  FORMED  AT  POINT  7 
LOAD  FACTOR  ADDITIONAL  CUMULATIVE 


STAGE I  S>  6.8895929E  01  2.663013AE  03 


DEFLECTION 

ADDITIONAL 

cumulative 

Xl 

1) 

-6.312lO2Ae-04 

-l  * 1415524E-02 

X< 

2) 

-3. 760401 5C-04 

-8.2191770E-03 

X( 

3) 

-6.3121024E-04 

•1  •  1415524C-02 

x< 

A) 

-5  *  372001 5E“0A 

-9.3607293E-03 

XI 

5) 

2.686000ef ”0A 

-3.652960OE-O3 

x( 

6) 

-5.3720015E-0A 

-9.3607296E-03 

xi 

7) 

9.7193713E-03 

1 .01217641-01 

X( 

6) 

1 .67^274  IE-02 

2.3896497E-01 

MO 

*£NT 

additional 

CUMULATIVE 

PLA5  MOM 

Ml 

1) 

»  i  > .  r.  a  • 

5999.9999 

6*14.  0000 

M( 

2) 

183.7225 

2301.3696 

6*04.4000 

M( 

3) 

111.2005 

3287.6711 

6000. 0000 

M  1 

A) 

145.0441 

4027.3973 

6*00. noon 

Ml 

5) 

-  r  .  r  ,  " 

5999.9999 

64C4.0000 

Ml 

6) 

321.5144 

4027.3969 

6ftoo. onon 

Ml 

7) 

176.4703 

5999.9999 

6400.0000 

Ml 

6) 

*  • 

. 

5999.9999 

6404.0000 

Ml 

91 

*  *•>  »  1  “  V 

5999.9999 

6404.4000 

Ml 

10) 

183.7225 

2301.3696 

641)4,  «oon 

Ml 

11) 

111.2005 

3267.6712 

6404.0000 

Ml 

12) 

145.0441 

4027.3971 

6400.4000 

Ml 

13) 

•294.9229 

-5569.0405 

6404.4000 

Ml 

U) 

•248.9923 

-5013.6979 

6444,0000 

Ml 

15) 

-248.9923 

•5013.6979 

6400.0000 

Ml 

16) 

-294.9229 

-5569.0404 

6404,4000 

Ml 

IT) 

•145.044 

•4027.3968 

6004.4000 

Ml 

16> 

•  o ,  (" 

-2999.999? 

6400,400m 

Ml 

19) 

-2999.9996 

6404.41)00 

Ml 

20) 

-145.0440 

-4027.3969 

6044.4(10*) 

hinge  ROTATIONS 


AT  POINT! 

.m  .  k 

-1.19 15306E-03 

AT  POINT ( 

£  . 

.«'-**  '  E-39 

AT  POINT! 

j! 

•  f  E-39 

AT  POINTI 

.  E-39 

AT  POINT! 

5) 

-2.T3972OC-03 

AT  POINT! 

r.  ;r,  .  E-19 

AT  POINT! 

V 

•  Mi"  E-39 

AT  POINT! 

6) 

.4.5662099E-OS 

AT  POINT! 

9’ 

-l*lAl550|E-03 

A?  POINT' 

10) 

.  .  t-39 

AT  POINT! 

v.E-39 

AT  POINT! 

(-39 

AT  POJNIl 

•  .  >v  ^  •  1-39 

AT  POINT! 

U' 

■ ,«  ,w ...  •  e-39 

AT  POINT! 

•  .  \.r*  •  ’,1-39 

AT  POINT! 

u» 

•■•pci.?*-  1-39 

AT  POINT! 

IT) 

1  .‘MV*'  1-19 

AT  POINT! 

U> 

...C"0C4',vI-J9 

AT  POINT! 

If) 

AT  POINTI 

20) 

•••it- of  44.  (.39 

-'W 
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PLASTIC  HINGE  No.  6  FORKED  At  POINT  13 

load  factor  additional'  cumulative 


2.7574467E  03 

cumulative 

-1 .241134BE-.02 
-8.5106377E-03 
-1 .2411348E-02 
-1  *0283688002 
-3.1914886E-03 
-1 .0283688E-02 
1.1465721E-01 
2.6418A42E-01 


cumulative 

5999.9999 
2489.3617 
3510.6385 
4276.5962 

5999.9999 
4595.7452 

5999.9999 

5999.9999 

5999.9999 
2489.3617 
3510. 63P6 
4276.5959 
-5999.9999 
-5297.8721 
-5297.8721 
-5999.9999 
-4276.5956 
-2999.9997 
-2999.9998 
-4276.5957 


HINGE  ROTATIONS 
AT  point!  1)  -2.6595752E-03 
AT  POINT!  2)  o.oCUOnnoE-39 
AT  POINT!  3)  O.OCOO.-,f:CE-39 
AT  POINT!  4)  ?  •oOC'l-CoE-39 

AT  POINT!  5)  -4 .6099302E-03 
AT  POINT!  6'  ;  *  j  0  0  0  •’  C  v  t  -  3  9 

AT  POINT!  7'  -8 . 8652689E-04 
AT  POINT!  -6 ,2056759E“°3 

AT  POINT!  9 I  -2. 6595752E-03 
AT  POINT*  10'  r.OOOO^CE-39 
AT  POINT!  11)  r.'.0CG0"^t.E-39 

AT  POINT!  12'  •• .  .lOO0r,,.UE“39 

AT  POINT!  13)  '■j.0C'?0r>r")E-39 

AT  POINT!  14)  ..  .’iGOO“:.';!E-39 

AT  POINT!  15'  f>  •jOOO*(',.-E-39 
AT  POINT!  16'  o.'ior,.'^or.E-39 

AT  POINT!  17)  >j.OOGG',OuE-39 
AT  POINT!  18>  ",0000nni>E-39 

AT  POINT!  19'  S..JCCC0'?iC-39 

AT  POINT!  20'  -j.00i.or  •.■•:,E-39 


PLAS  MOM 

6OUO.O0OO 
6000.0000 
6000. 0000 
6000.0000 
6000.0000 
6000,0000 
6000.0000 
6noo.ooon 
6000. ooco 
6000 .0000 
6000.0000 
6000.0000 
6000.0000 
6000 ,0000 
6000.0000 
6000.0000 
60.00,0000 
6000. 000 
6000.000 
6000  .noon 
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D 


plastic  HitfGE  mo*  7  formed  point  u 
load  FACTOR  ADDITIONAL . CllMU'LATiVr 


STAGE!  7)  1.9151234E-05 

2.7574467E  03 

deflecti0n  additional 

cumulative 

X  ( 

1) 

-3*  7942&77e**I0 

-1.24U348E-0? 

X! 

2) 

-1.3346932E-1Q 

-P.5l06.378F  -f)3 

X  ( 

3) 

-2.7234954E-10 

-l.P'.HlAAf-OP 

X  ( 

4) 

-2*  18904«5E'-10 

-1 «0?83688F-0? 

x< 

5) 

1  .U37S70E-10 

-3.  i9i<*e«6r-o3 

X! 

6) 

-2.54S9793E-10 

-1 .02836891-02 

X  ( 

7) 

3*7996672£-09 

1 1 14647??r-0l 

X  ( 

8) 

7 • o5fe 1 8 56E-09 

?.841844?f.-01 

MOMENT 

additional 

cumulative 

*LAS  MOM 

M  ( 

1) 

-f'.OMH, 

5999.9999 

6609.0000 

M! 

2) 

0  ,  o  r  0  r. 

2489.3617 

6000.0000 

M( 

3) 

-O.'Wv 

3510.6385 

6000. 6000 

M< 

4) 

C)  *  'T  i' ■ 1 

4276.5962 

6009.000(1 

M  ( 

5) 

-0,0'- -Jv-- 

5999.9999 

6000.9000 

Ml 

6) 

C.o.-u 

4595.7453 

6000. 6000 

M  ( 

71 

0.0 

5999.9999 

6000. 60Q0 

M  ( 

8) 

^  o.  -. 

5999.9999 

6OQ0. 9000 

M  ( 

9) 

-0  .  '} ! ' 

5999.9999 

6009.0000 

Ml 

10) 

o.o*  ■:>  l 

2489.3618 

6009.9000 

Ml 

11) 

o  .  0  r  n  1 

3510.6386 

6000.9000 

Ml 

12) 

4276.5960 

6500.0000 

Ml 

13) 

-o.o-.n: 

-5999.9999 

6onotooon 

Ml 

14) 

— l  .  'j  n  ’  i 

-5297,8721 

6000.0000 

Ml 

15) 

-n.o-.oi 

-5297.8721 

6900.0000 

Ml 

16) 

-n  .0~"1 

-5999.9999 

6600.9000 

Ml 

17) 

-0.0: 

-4276.5956 

6906.9000 

Ml 

18) 

0 . 0  r  j  i; 

-2999*9997 

6oon. ooon 

Ml 

19) 

—0 .  9  ^  *  <  o 

-2999.9998 

6606,0000 

Ml 

20) 

-O.0r.Cl 

-4276.5958 

6900. O000 

HINGE  ROTATIONS 

AT 

POINT ( 

1)  -2.6595755E-03 

AT 

POINT ( 

2 1  .!.iaC0"0ve-39 

AT 

POINT ( 

3 1  r .  ♦  U  C«  0  0  C  0  *  ■  I  -  3  9 

AT 

POINTI 

4>  o,  0C(.-0ni;!iE-39 

at 

POINT  ( 

5  >  -4 . 60 99306 E -D 3 

AT 

POINT ( 

6>  i.  dOiJOr  urE-39 

AT 

point'- 

7)  -8.8652689E-04 

at 

PO  INTI 

8)  -6 .2056761 E-03 

AT 

POINT ( 

9>  -2. 6595755E-03 

AT 

POINT  f 

10)  '•  •  'jOQO  *  tK'E-39 

AT 

POINT  ( 

ID  'i  •00000fi|<  E  -39 

AT 

POINT! 

12)  u*  0(-0!'»p«;:>E-39 

AT 

POINT' 

1 3 1  '■ .  0  c  0  0  i>  < :  E  -  3  9 

AT 

POINT' 

1  4  •  1 ' .  u  <"•  1  ■  0  IT  0  U  E  -  3  9 

AT 

POINT! 

15-’  ”  •  0  0  u 9  A  C  E  *  3  9 

AT 

POINT ( 

1 6  )  •  (i  0  \j  ()  *  i:  '1 £  -  3  9 

AT 

POINT  ! 

17 »  ■..cooo.70uE-39 

AT 

POINT ( 

lfl'  I..0O007IJL-X-39 

AT 

POINT  f 

19!  '.’.OCofno:E-39 

AT 

POINT ( 

20)  •».00«)f:'*.0CE-39 
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PLASTIC  HINGE  NO.  8  FORMED  At  POINT  6 


LOAD  FACTOR  additional 
STAGE  I  8)  1.4781631E  02 


DEFLECTION  ADDITIONAL 


X  ( 

1) 

-5-  13250.95E-03 

x< 

2) 

-2.6004716E-03 

X  ( 

3) 

-5 •  1 325095E-0  3 

X  ( 

4) 

-2.8742050E-03 

X  ( 

5) 

1.4371025E-03 

X  ( 

6) 

-2.S742050E-03 

X  ( 

7) 

5-2009434E-02 

XI 

8) 

9.5464680E-O2 

MOMENT 

additional 

M( 

1) 

.  0  ,  ,'i - •_ 

M< 

2) 

36.9543 

M( 

3) 

-36.9541 

M< 

4) 

776.0355 

M( 

5) 

-  o.c 

M  ( 

6) 

1404.2547 

M< 

7) 

C  .  0  °  1 

Ml 

8) 

0.0*  Ilf 

M  ( 

9) 

- 1  •  .  0  9  *  i. 

M< 

10) 

36.9543 

Ml 

11  ) 

-36.9541 

M  I 

1?) 

776.0355 

Ml 

13) 

*•  v.  •  1 !  *  «.  i  •  » 

M( 

14) 

-702.1273 

Ml 

15) 

-702.1273 

Ml 

16) 

-  -it*.' 

M( 

17) 

-776.0353 

M  ( 

18) 

r  .  ■)  f  r>  r 

m( 

19) 

'•.7- 

Ml 

20) 

-776.0353 

HINGE  ROTATIONS 

AT 

POINT' 

1) 

-7.8947357E-03 

AT 

POINT  < 

2) 

0 .  'jC CO nO.'.£ -39 

AT 

POINT ( 

3> 

*  i  *  U  D  0  C  ,  't  ■  )  i  C  —  J  9 

AT 

POINT! 

4 ) 

n.ononn:j.;E-39 

AT 

POINT ( 

5) 

-1  • 1 1 1 1 1 1 0E-02 

AT 

POINT ( 

6) 

r*.  )0 c n 1> ;.-E - 3 9 

AT 

POINT ( 

7) 

-2. 63158006-03 

AT 

POINT ( 

85 

-1  .  1988303E-02 

AT 

POINT ( 

9> 

-  7  •  894  7357E-0  3 

AT 

POINT ( 

10) 

. .  u  0  *  n  j  v.  E  -  3  9 

AT 

POINT ( 

IP 

*  •  uOuG*  00E-39 

AT' 

POINT ! 

12> 

« < .  j  0  0  0  *  0  u  l  -  3  9 

AT 

point' 

13) 

6.4327457E-03 

AT 

po  in  r ( 

14) 

i  E -39 

AT 

POINT ( 

15» 

0OuO'V.i  E“39 

AT 

POINT i 

1  6 ) 

6.43274566-03 

AT 

POINT ( 

17) 

*;.0000n-K*E-39 

AT 

point  ( 

18) 

•  t .  ijOuOO'i'.'E-  39 

AT 

POINT ( 

19) 

o.nODG*  tce-39 

AT 

POINT ( 

20) 

*.0'VjC*HE-39 

cumulative 

2.9052630E  03 

CUMULATIVE 
-1.7$4385?F-02 
-1.1UU09E-02 
-1.7543857E-02 
-1 .3157893E-02 
-1.7543859E-03 
-1 j  3157894E-02 
1 ;6666665E-0l 
3t59649l0E-01 


CUMULATIVE 

PLAS  MOM 

5999.9999 

6*o*.oooc 

2526.3160 

6*0* .0000 

3473.6844 

600*. 0000 

5052.6317 

60oc .oonn 

5999.9999 

6000 , 0000 

5999.9999 

6*00.0000 

5999.9999 

6000,0000 

5999.9999 

6*00.0000 

5999.9999 

6*0*.  0()00 

2526.3160 

6*0*. OQOO 

3473.6845 

6*00. *000 

5052.6315 

6000.0000 

-5999.9999 

6*0*. 0000 

-5999.9994 

6*00.000* 

-5999.9995 

600*. 0000 

-5999.9999 

6*00,0000 

-5052.6309 

6*0*. *000 

-2999.9996 

6*co.*ono 

-2999.9998 

6*0* , 0000 

-5052.6311 

6*0°  .0000 
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PLASTIC  HINGE  No.  9  FORMED  At  POINT  4 
LOAD  FACTOR  ADDITIONAL . cumulative 


STAGE ?  9) 

9.4736889E  01 

DEFLECTION 

ADDITIONAL 

X  ( 

1) 

-7*456i370E“03  .. 

X( 

2) 

-O.GOCO*1'  E-39 

X  ( 

3) 

-7.4561370E-03 

XI 

4) 

-3 .506770 1 E-03 

X( 

5) 

1 . 7543851 E-03 

X( 

6) 

-3* 50 87701 £“03 

X( 

7) 

8.3333299E-0* 

X  1 

8) 

1.4035082 £-01 

MOMENT 

ADDITIONAL 

M{ 

1) 

-f.QroO 

Mt 

2) 

473.6841 

M< 

3) 

-473*6837 

M  ( 

4) 

9^7.3683 

M< 

5) 

—  0  «  Oi  ’i'.'' 

M  ( 

6) 

_  r .  n  *>  i.i  j 

Ml 

7) 

'1.0  '■  ') 

M< 

8) 

f-.DOliU 

M< 

9) 

-C  .  OrO'.i 

M  ( 

10) 

473*6841 

M  ( 

11) 

-473.6837 

M  ( 

12) 

947.3fc83 

Ml 

13) 

—  0  •  0  0  0  0 

Ml 

14) 

•O.CrOc 

Ml 

15) 

«Q.0**0 

Ml 

18) 

-  0 . *  •'  c  0 

Ml 

17) 

-947.3679 

Ml 

18) 

0.0000 

M  < 

19) 

0  *  3 n<).; 

M  < 

20) 

-947.3679 

HINGE  ROTATIONS 

AT 

POINT' 

1 ) 

-1 .6666662E-02 

AT 

POINT < 

2) 

,j .nOOCnOC  E-39 

AT 

POINT  I 

3> 

-  .  ^  • 1  j  C  —  3  9 

AT 

POINT? 

4) 

-*  ->E-39 

AT 

POINT ( 

5) 

-1 , 944443  9E-02 

AT 

point? 

6’ 

-8* 333330CE-03 

AT 

point? 

7) 

-8  *333331 9E-03 

AT 

POINT? 

8) 

-1 . 944  4  440E -02 

AT 

POINT ( 

9> 

-1 • 6666662E-02 

at 

point? 

10) 

.  !.•  0  0  0  <•  1  •  1 E  -  3  9 

at 

POINT? 

ID 

fi.0nu0P'JyE-39 

AT 

POINT? 

12> 

r .  0  7  0  0  0  •:)  r.  E-39 

at 

point? 

13) 

1 .3888883E-02 

AT 

POINT ( 

14) 

.*1 . 0 *  C  0  *  *  G  E  -  3  9 

AT 

POINT  ? 

15) 

' .  0  0  0 ' )  ■  *  0  £  -  3  9 

AT 

POINT ( 

16) 

1.3888882E-02 

AT 

POINT? 

17) 

0.0'’CO'' j-:'E-39 

AT 

POINT? 

IV 

■1.0  0  n  0  <-»  0  L*  E  -  3  9 

AT 

POINT? 

19) 

7  •  0  0  0  0  0  0  U  £  -  3  9 

AT 

POINT ( 

20) 

O9u0r  70E-39 

2.9999998E  03 

CUMULATIVE 

-2,A999994r-0i 

-i • 1 11 1 109E-02 

-2.4999994E-02 
-I.6666663E-Q2 
-8,440] lOgE-10 
-1 .666666AE-02 
2.4999995E-Q1 
4-9999992E-01 


CUMULATIVE 

PLAS  MOM 

5999.9999 

6000. 0000 

3000 , 0001 

600,'.ooon 

3000.0007 

600*. 0000 

5999.9999 

6000. oonn 

5999.9999 

6000.0000 

5999.9999 

6000.0000 

5999.9999 

6000. oonn 

5999.9999 

6000,0000 

5999.9999 

6000,0000 

3000. COOl 

600^.0000 

3(j00 . 0?)08 

6000.0000 

5999.9998 

6000.0000 

-5999.9999 

6000.0000 

-5999.9994 

6000. oono 

-5999.9995 

6000.0000 

-5999.9999 

6000. 000O 

-5999.9988 

6*00,0000 

-2999.9996 

6*00.0000 

-2999.999 $ 

6000,0000 

-5999.9990 

6*00.0000 
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PLASTIC  HINGE  HO.  10  FORMED  AT  POINT  12 

ioao  factor  additional'  cumulative 


STAGE (  10) 

1 » 2207fl  35E-05 

2 » 9999998E  03 

DEFLECTION 

additional 

CUMULATIVE 

X< 

1) 

-1.AA11O7.AE-09 

-2. A999996E-02 

X( 

2) 

■T  ■  0  0  0  0 ' 0  E  -  3  9 

-1.1I1U09E-02 

x< 

3) 

-1.55A13S3E-09 

-2.A999995E-02 

X  ( 

A) 

-1.1302803E-10 

-1 .6666663E-02 

X  t 

5) 

2.2605605E-10 

-6. 1795503E-10 

X  ( 

6) 

-7.9U9615E-10 

-1 , 6666664F-02 

X  I 

7) 

i  .610649AE-08 

2.A999996E-01 

x( 

8) 

2.88221A8E-08 

A.9999994E-01 

MOMFNT 

ADDITIONAL 

CUMULATIVE 

PLAS  MOM 

MI 

1) 

-  n  •  0  " 

5999.9999 

6000,0000 

M  ( 

2) 

0 . 0  " 1 

3000.0002 

6000,9000 

M  ( 

3) 

-o.oooi 

3000.0006 

6000. *5000 

Ml 

A) 

0, 

5999.9999 

6000.0000 

M( 

5) 

—  C  •  9  "O'.. 

5999.9999 

60(10.0*100 

M  ( 

6) 

•0.0  “  ;l 1 1 

5999.9999 

6000.0000 

M( 

7) 

O.rji'o  . 

5999.9999 

6000.0000 

Ml 

8) 

0  •  0  0  sj  o 

5999.9999 

6000,0090 

M< 

0) 

-  0 . w 

5999.9999 

6000.0000 

M  I 

10) 

u  .  0  -  • 

300^,0002 

6000. OOOO 

Ml 

ID 

-  ?  • .  0  0  0 

3000.0008 

6nc^.0000 

Ml 

12) 

0.00  0  2 

5999.9999 

6nno.000n 

Ml 

13) 

•  0.0**  0  v,1 

-5999.9999 

6000.0000 

Ml 

1 A ) 

•o.n-m.i 

-5999. 999A 

6000,0000 

Ml 

15) 

-O.OoO'.- 

-5999.9995 

6000,0000 

Ml 

16) 

*  0  •  u  t  j  1  u 

-5999.9999 

6000. 0000 

Ml 

17) 

-0  ,  Or r  ; 

-5999.9988 

6000, «oon 

Ml 

18) 

O.O^ni 

-2999.9995 

60(10.0000 

Ml 

19) 

-  o .  o  r  o  1 

-2999.9998 

6000.^000 

Ml 

20) 

•  0 .0  '*•  0  2 

-5999.9992 

6000.0000 

HINGE  ROTATIONS 

AT 

POINT  ( 

I ) 

-1 .6666663E-02 

AT 

POINT ( 

2) 

IJ.0000O00E-39 

AT 

POINT  * 

3' 

?  0  o  i)  n  -i  0  E  “  3  9 

AT 

POINT* 

A) 

•,  O';  CO. VOCE -3  9 

AT 

POINT ( 

5) 

-1 .9A4A4Alf-02 

AT 

POINT ( 

6> 

•  8*  333331  AE-03 

AT 

POINT l 

7) 

-8.3333333E-03 

AT 

point ( 

8) 

-1 .9A<*ALA2E-02 

AT 

POINT ( 

9’ 

-D6666663E-02 

AT 

POINT* 

10) 

c  .ooro*1''-  E-39 

AT 

POINT ( 

ID 

O.OOii'^f'LE-39 

AT 

POINT  * 

1  2 ) 

5  •  0 CuO*1 ;.! i  E - 3 9 

AT 

POINT* 

13) 

1 .388888AE-02 

AT 

POINT  * 

1A> 

<- ,U00i)-'  TuE-39 

AT 

POINT ( 

15) 

■  j .  iT  0  c  0 i*'  ' '  E  —  3  9 

AT 

POINT* 

U) 

1.388888AE-02 

AT 

POINT ( 

17) 

r. ,  000ri:’>r!,'E-39 

AT 

POINT ( 

18) 

c  *  C'ft00r  j uE-39 

AT 

POINT* 

19) 

0  •  OOCUfM'C  E-39 

AT 

POINT* 

20) 

OOO'T '•OQE-39 
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APPENDIX  C 


COMPUTER  INPUT  AND  OUTPUT  LISTINGS 
FOR  SAMPLE  PROBLEM  1 


Preceding  Page  Blank 
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DISTANCES  FROM  FIRST  HALL  5C • * 15C. * 200.  FEET 
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